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SPECTRAL INVARIANT OF THE ZETA FUNCTION 
OF THE LAPLACIAN ON S*r-} 


N. STHANUMOORTHY 


Ramanujan Institute for Advanced Study in Mathematics 
University of Madras, Madras 600 005 


(Dedicated to Professor T. S. Bhanu Murthy on his 60th Birthday) 
(Received 8 June 1987; after revision 26 October 1987) 


The aim of this paper is to compute a spectral invariant of the zeta function 
0 (4, s) at s = Oof the Laplace Beltrami operation A acting on forms of 
degree 2 on (4r — 1) dimensional sphere S‘7~!, 


1. INTRODUCTION 
Let A be an elliptic pseudo-differential operator of positive order m on acom- 
pact n-dimensional manifold YX. If A is self adjoint and positive, its eigen values are 
A > 0. We define its Zeta function by 


€¢ (A, s) = Trace A~* = X= A-* (Atiyah ef al), 
»A>0 


Here each eigen value is repeated as many times as its multiplicity. This series con- 
verges for Re (s) > es (Atiyah et al.*) and hence gives a holomorphic function of 
the complex variable s in this half plane. Also ¢ (A, s) can be analytically continued to 
the whole s-plane as a meromorphic function with simple poles. On the other hand, 
for any elliptic self adjoint pseudo differential operator B of positive order m ona 
compact manifold X, we define its eta-function as 


n(B,s) = = Signala|~*. 
40 


Here also A runs over the eigen values of B and each eigen value repeats as many 
times as its multiplicity. But here as the operator B is not positive, its eigen values 
may be positive or negative and hence each eigen value is taken with its sign. When 
the operator B is also positive, we have 


¢ (B, s) = 1(B, 5). 


The real valued invariants of the metric satisfying the condition that it is a con- 
tinuous function of the metric can be obtained by evaluating ¢ at some point s where 
it is known to be finite. For positive self adjoint elliptic operators, the zeta function 
defined above have finite values at s = O by the results in Seeley®. Atiyah ef al.*, gave 
two methods to prove the finiteness of » at s = 0 for the special operator B,*’. Using 
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Cobordism theory one can prove the finiteness of 4 (0) and by using invariant theory, 
1 (0) can be computed explicitly for B,*’. But in this paper we show a general method 
of analytic continuation to compute ¢ (A, s) ats = 0 for the Laplacian 4 and compute 
it for A acting on 2-forms on S*’-'. The importance of this paper lies in the analytic 


co 
continuation that we have used for the function of the form 2g (n) {f (n)}-’, where 


f (n) and g (n) are two polynomials of finite degrees. The method of analytic continu- 
tion that we have used in this paper can be used for all similar series. 


2. EIGENVALUES OF A ON 2-FORMS ON S4'~? 


When G is a compact connected Lie group and Kis a closed subgroup, we con- 
sider the quotient space M = G/K. The Laplace Beltrami operator or Laplacian is 
A = d83 + 3d, where d is exterior differentiation and Sis operator adjoint tod. Aisa 
self adjoint and elliptic differential operator on C® (A?M) for each p, Here C®(A?M) 
is the vector space of smooth sections of the vector bundle A?M, the pth exterior 
power of the complexified cotangent bundle of M. The set of eigen values of A on 
C~” (AM) is a discrete set of real numbers. 


When G is a compact semisimple Lie group and B is the killing form of the Lie 
algebra CG of G, the Casimir element isC = cl X,.X, where {Xj,....Xn} is a 
151, jSN 
basis of Zand c’ = (B (X;, X,))-1. When (G, K) is a compact symmetric pair with a 
compact connected semisimple Lie group G, we have the Cartan decomposition 


SF=kOn, 


where k is the Lie algebra of K and m is the orthogonal complement to k in TZ with 
respect to the Killing form. Restricting the Killing form sign changed to m, we get a 
G-invariant Riemannian metric on M = G/K andA = —C (Ikeda and Taniguchi?). So 
using Proposition 2.2 of Ikeda and Taniguchi’, the eigenvalues of A can be written for 
S" = SO (n+1)/SO (n). We are interested in the Laplacian A acting on 2-forms on 
S41, Let t be a Cartan subaigebra of Fand Ai, ..,Asr be linear forms ont. Any 
dominant integral form A in G = SO (4r) with respect to t is uniquely expresséd as 


A => Ky Ai + eee + Ker Aor, 
where k,,...,ke, are integers such that 


hy PK a | Re |. 


Then using the theorem 4.2 of Ikeda and Tani ; 
: uchi*® the foll 
easily stated, é Owing theorem can be 


Theorem 1—If p == 2 and Ap is the hi 


. ghest weight of their i : 
P intervening in C? (A2S*"-2) then 8 reducible representation 


Ap = (kK +1) A, +A, or(k+1) a1 + A, + A; for r > 2. 
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Moreover with respect to CES) times the Killing form, the eigen values of A are 
given by 
(Kk + 2)(k + 4r — 2) for Ap = (kK + 1) A, +A, 


and 
(kK + 3)(k + 4r —3) for Ap = (K + 1) A, + Ao + As. 


Here k runs over all non-negative integers. Moreover the multiplicity pe of the 
above in C® (A? S*"~") is exactly one in the above two cases. 


The multiplicity of an eigen value is the product of the multiplicity ue and the 
dimension of the representation P. The dimension of the representation P can be com- 
puted using Weyl’s formula : 


a < Ap + 3,a> 
>t : 


i Oe: 
a>0 


Here « > 0 are the positive roots in SO (4r) and 5 is half the sum of the positive roots 
in SO (4r). 


Hence we get the following table : 


Ae Eigen value 
(A+ 1)A, + A, (kK + 2)(k + 4r — 2) 
(A+ 1)A, +A, +A; (k + 3) (kK + 4r — 3) 


Multiplicity 


2(k + 1) (kK + 4r — 1)(k 4+ 2r) Gamera 
4 -3)h+2) 0° 4r — 4 

(k + 1) (kK + 2r) (k + 4r — 1)(4r — 3) Gee) 
Bee hk 43). 4r — 3 


Here k runs over all non-negative integers. 


3, ZeTA FUNCTION OF THE LAPLACIAN A FOR S‘*~? 


As we consider the positive operator A acting on 2-forms on S“-1, we get 
n (A, s) = ¢ (A, 5). 


2 (n +1) (n + 4r — I) (n +: 2r) (n + 4r — 3 
“OD = SS Sea eel ar ea 
x {(n + 2) (n + 4r — 2) }°° 


(equation continued on p. 410) 
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(2 + 1) (n # 2r) (a -Padr — 1).4r 3) 
" (n + 3)(n + 4r — 3) 


n=0 


(Ets?) e mete af 


we have to compute ¢ (A, s) at s = 0 by a method of analytic continuation. We denote 
this value by ¢ (A, 0). 


We explain below the details of the analytic continuation which we use to compute the 
spectral invariant of the zeta function for S*’~*. 


Let 
a = 8 {fy 


Here g (n) and f (”) are two primitive polynomials of finite degrees in 7. Let us assume 
that the degree of f (7) be k. Then this series S has analytic continuation as 


el g (n) { (n*® + P (n))-*§ —n-*s} + = g (n) n-* 
a= ni 


in the entire complex plane. 


Here f(n) = n* + P (n), P (n) being a polynomial of degree k —1. 








Let 
(I) = = g(n){ (nk + P(n))-? —n-*} 
n<c 
ier rele (eee ape || 
(III) = z g (n) n-*s, 


Here the positive number C is chosen such that 


oss 





~<a 


Then S has analytic continuation in the entire complex plane with 


S = (I) + (Il) + (IID. 


(1) is an entire function whose value at 5 = O%s.00 ind); C= 3 By * can be ex- 
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panded using binomial theorem because | Po < 1. Hence 
(1) = SS gina | —s he += eae? ( = — ... to cof, 


n>C 


As s >0Owe haves¢(s + 1) > 1, © being the ordinary Riemann zeta function’. 
So when S — 0, (II) gives some constants due to first few terms and all the other 
terms in (II) will tend to zero. Morever the sum in (II) taken over any finite rectangle 
tends to zero as s > 0. Let g (n) be a polynomial of degree ¢ such that 


q 
g(n) = & an‘ where a (for 0 <i < q) 
i=O 


are constants with a, = |. 


Then 
(I) = © at (ks — i). 
{-0 


So when s > 0, (III) will contribute some constants to ¢ (A, 0) and finally we get 


t (A, 9). 
We now compute the value ¢ (A, 0) for S*~?. 
Let 
2a ee (ee 2 dh) 
(4r — 3)(n + 2) 4r —4 (4r — 3)! °’ 


et . . 
where g, (n) = = an‘ with ay. = 1. Similarly 
iO 


we assume that 


(4r — 3)(n + 1)(n + 2r)(n + 4r cee) fi + 4r — *) 
(n + 3)(n + 4r — 3) 4r — 3 


] 


4r-—2 
where g2(n) = = bn‘ with by-. = 1. 
i~0o 


We also assume that 


fin) = (n+ 2) (0 + 4r — 2) 
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and 

f, (n).= (n + 2) (n + 4 — 3). 
So 

C (A, s) = 2r (4r — I + (ar = 1) (4r — 2) 
(at s°=— 0) 


+ (pees s) (L@"t 


l _ =g Pa 
+ gem Ds gs (0) (fa (0) (ats = 0) 
We first find the contribution to ¢ (A, 0) from 
8 = Fam(Amy 


This series has analytic continuation in the entire complex plane as 


S, = (I) + (I) + (IID) where 


li 


(hes Ss 81 (n) | (n? + 4rn + 8r — 4)-* —n-2s } 


n=l 


(I) = S gy (n)n-2s {a i wae 4 =] } and 


s=C+1 


(II) = SS 81 (n) n-*s, 


m=) 


Here C is chosen such that Ma BaRlint ig 
n* ‘ 


Remarks : When r = 


2 (S"), C can be chosen to be greater or equal to 9. W 
r = 3(S!'), C can be chosen ; . rs 


to be greater or equal to 13 and so on. 


(1) is an entire function whose value ats = Ois 0. 


co 
(II) = SS g, (n) nn ig AO Brie a 
n=C+1 n* 


(equation continued on p. 413) 
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- 2 
ee) al 


n2 
Let P; (n) = 4rn + 8r — 4 and 


4r+4f-—2 
£1 (n)(P, (n)) = eS Babel ear Sd. 


i=0 


Denoting by (II);, the ith term of (II), we get 


co 
(II), = —s 2 & (n) Pi (n) n-@st?) 


4r 
(I; (at s = 0) — —<¢ E aon E(2et.2 —ddars t= 0) 


= — $ dQ). 
Similarly 
(II), (at s = 0) = 3 ag,3), 
(IDs (at s = 0) = —% a,5) and so on. 
Finally 


(— 1)" 


(11) sr-1 (at a 0) = 2 (4r — 1) Q(4r—15 gr—3) 


and all other terms in (IJ) will tend to 0. Now 


io) 
(I) = = gi (n)n-* 
n=1 


4r—-2 
= x Qj C (2s — i) 
i=(0 





iy, a, CET a3B, 
(III) (at s = 0) = D 5 ate 
(=I) ers Bers 
. + 4r o 2 . 


Here &,, B.,...,B,r-, are Bernoulli’s numbers’. 


Similarly to find the contribution to ¢ (A, 0) from 


S: = = 82 (n) (fo (n))-’~ 


414 N, STHANUMOORTHY 


we assume that 


P,(n) = 4rn + 12r — 9 


and 
4r+t—2 


es (m) (Pp (2)! = "Eben mt for < tar = 1, 


So using the notations that we have used, the following theorem is completely 
established. 

Theorem 2—A spectral invariant of the zeta function ¢ (A, s) at s = 0 of the 
Laplace Beltrami operator A acting on 2-forms on S4"-! (r > 2) is 


r- or =) 
2 — (=I)! G(t9¢-1) ut >; (=I)! Gy-) Be \ 
(4r — 3)! DS 2t 2t - 
me | 


t=) 


ar-1 ah 
l ' (=1)! Detyge—1) ws > (= 1)! bur-) Be oa >| 
1 (4a — 4)1 reer 21 2 


t=] i=) 
wi a (16r? —1). 


Remarks; The author is thankful to the referee for pointing out an error in 
the eigen values and making some encouraging comments. Using the results of 
Sthanumoorthy*”’ some results on spectral invariants will be published shortly. 


ACKNOWLEDGEMENT 


The author thanks Professor K. Ramachandra of Tata Institute of Fundamental 
Research, Bombay for very useful discussions he had with him. From these discussions 
it was understood that Professor K. Mahler had long back obtained more general and 
deeper results on the convergence of such series and integrals. The author is also 
thankful to Professor S. Raghavan of Tata Institute of Fundamental Research, Bombay 
for some useful discussions. 


- 


REFERENCES 
1. M.F. Atiyah, V. K. Patodi, and I. M. Singer, Bull. Lond. Math. Soc. 5 (1973), 229-34. 
2. M. F. Atiyah, V. K. Patodi, and I. M. Singer, Math. Proc. Camb. Phil. Soc, 79 (1976), 71-99. 
3. A. Ikeda, and Y. Taniguchi, Osaka J. Math. 15 (1978), 515-46. 
4. K. Mahler, Mathematische Annalen, 100 Band (1928), 384-98, 
5. R.T. Seeley, Complex powers of an elliptic operator. Proc. Symposium in Pure Math. Vol. 10 
Amer. Math. Soc., 1967, pp. 288-307. 
6. N. Sthanumoorthy, Bull. Se. Math., 2¢ Serie, 108 (1984), 297-320, 
i N. Sthanumoorthy, Bull. Sc. Math., 2° Serie, 111 (1987), 201-27. 


E. C. Titchmarsh, The T, heory of the Riemann Zeta Function. Clarendon Press, Oxford, 1967 


Indian J. pure appl. Math., 19 (5) : 415-422, May 1988 


INCOMPLETE BLOCK DESIGNS OBTAINED BY THE GENERALISED 
ROW-JUXTAPOSITION AND GENERALISED COLUMN- 
CONCATENATION OF INCIDENCE MATRICES 


G. A. PATWARDHAN 
AND 


SHAILAJA SHARMA 
Department of Mathematics, Indian Institute of Technology, Bombay 400 076 


(Received 15 January 1987; after revision 12 October 1987) 


We define the concepts of ‘generalized row-juxtaposition’ and ‘generalized 
column-concatenation’ of two matrices. By taking these matrices as incidence 
matrices of block designs, several new constructions of binary and n-ary 
(n>2) incomplete block designs have been obtained. 


1. INTRODUCTION 


Taking the clue from Helleseth and Van Tilborg*, we have defined, in the 
following, the concepts of ‘generalized row-juxtaposition’ and ‘generalized column- 
concatenation’ of matrices. A reference to what we have termed ‘generalised column- 
concatenation’ has been made in Helleseth and Van Tilborg’. 


The generalized juxtaposition of the rows and concatenation of the columns of 
(0, 1)-matrices provide some new constructions of 2-PBIBDS, 3-PBIBDS, partially 
balanced n-ary designs and balanced block designs. Eleven such constructions have 
been given in section 3 through various theorems and corollaries. 


2. DEFINITIONS 


Let A and B be arbitrary matrices of orders m X n and p xX 4 respectively, over 
a field F. 


Let 
i Bie ak 
and 
aj; l1Sjgn 
denote respectively the row and column vectors of A. Let 


So he iP 
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and 
bj l<j<q 
denote the corresponding vectors of B. 


Given the ordered pair (A, B) of matrices, the generalised row-juxtaposition 2 
the pair (A, B), denoted G (A | B), is a function that juxtaposes all possible ae 
rows, one each from A and B taken in that order, giving rise to a new matrix G (A | 


of order mp x (n+ q): 


Hi 
r, r, r, . Fy m ™m 
A|B)= 
G(A|B) Cae. ia 


Analogously, given the ordered pair (A, B), the generalized column-concatenation of 
(A, B), denoted G (A/B) concatenates all possible pairs of columns, one each from A 
and B taken in that order, giving rise to a new matrix G (A/B) of order (m + p) Xnq. 


Re ag 
bisenbe 


The following standard abbreviations have been used : 


ay Peat 


Ao... Ae 
IS eit ve 


b, sale 











SBIBD _: symmetric balanced incomplete block design 


n-PBIBD : n-associate class partially balanced incomplete block design 


BBD : balanced block design 
SRGDD_ : semi-regular group divisible design 
S-PBnD : S-associate class partially balanced n-ary design 
C,(A) : Hasse-Minkowski invariant of A 
A (v, b,r,k, A) : BIBD with parameters y, b, r,.k, A having A as its incidence 
matrix. 


3. MAIN RESULTS 


The following theorems can all be proved by the actual construction of the 


matrices and computation of the design parameters indicated. Most of the proofs are 
straight-forward, hence omitted. 


| Theorem 3.1\—G (A | B), where A (1, 51, 71, ky, a3) and B(v,, bo, re, k 
BIBDs, is a 3-PBIBD with rectangular association sche 
and the following parameters : 


9, A.) are 
me, two unequal block sizes 


ace "1% B= db) +b, R= +r, +7; 
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Ki = k, v2, K, = ke v,; 
A. = rr +A, A, = re +A, As = Ar + Az 
m = ve— 1, m2 =v, — 1, my = (vy, — 1) (v2 —1) 


there being 5, blocks of size K, and 6, blocks of size Ky. 


Corollary 3.1.1—In the above construction, when A = B, G(A | A) represents 
a 2-PBIBD with Z»-association scheme, with parameters 


eet pa Be 20) Ro -2r, Ki = vk: 
A, = r+a, A, = 2); 
mn = 2(v—1), nm =(v — 1)*. 


Raghavarao‘ (p. 159) has given a construction fora 2-PBIBD with the above para- 
meters, wherein he is required to use a BIBD and also an association scheme of Latin 
square type for constructing this 2-PBIBD. Our construction requires only a BIBD 
A (v, 6, r, k, 4) and nothing else. 


Corollary 3.1.2—Each submatrix of G(A | B) with constant column sum consti- 
tutes a 2-PBIBD of the singular group divisible type. The parameters of these designs 
are 

(V, Voy Dy, 71, Va Ky 3 71, Ay 3 Ve — 1, (v1 — 1) vy) 
and 

(04 Ve, Dy, Ta, Vi Ka; Pe, Ags %,'— 1, (, — 1) ¥,) 
respectively. 
Since 

G (AT | BT)? = G (A/B), 
we also have 


Corollary 3.1.3—The transpose of G (A/B) is a 3-PBIBD with rectangular associa- 
tion scheme and unequal block sizes, provided that A (»,, 5,, r,, ki, A) and B (v2, b,, 
r,, ky, Az) are linked block designs with respective block intersection numbers #1 and 


fy. Its parameters are 
V=b6,b,, B=\y,+%, Rh +h; Ki = bn, Ky = Or: 
A, = ki + pe, A, = ke + ph, As = pa + pos mm = 4b, — | 
n, = b, — 1, nz = (b, — 1) (2 — 1). 


Corollary 3.1.4—If A = B in Corollary 3.1.3, then a 2-PBIBD with L.-association 
scheme can always be constructed, whose parameters are V = b?, B = 2v, R = 2k, 
K = br; A,=>k+ 4, A= 2p; 2, = 2 (6 — 1 ng = (6b — 1}*. 
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BDS, then 
Theorem 3.2—If A (v1, bi, 71, ki, Ax) and B (Yo, bey Ie, Kas jap ee a Kea 
G (A/B) is a BBD with two different numbers of replicates as defined by 
Its parameters are 


V=y+,, B= b, by; Ry = br, Re = bi Fe 
K=k, + ky; Ay = 6, A;, As =b1 A, Aye = 11 Te. 
This real is corroborated by the structure of G(A/B). [G (A/B)}’, which is 
bcd At oeel 
i 


ye | 


ry re . se - 


b, . BBT 





Corollary 3.2.1.—When A = B in the above theorem, the BBD becomes a 
SRGDD with parameters 


V = 2v, B= D8, R= rb, K = 2k 
A, = ba, A,= r*3 nj = v — 1, n, = ¥. 


Corollary 3.2.2—The transpose of G (A | B) isa BBD with two different numbers 
of replicates, provided that A and B are both linked block designs with respective 
block intersection numbers ;; and py». Its parameters are 


V = b, + b,, B=v»,;3 Ri = vk; Ry = v; ko; K=r, ite 
A, = v2 1, A, = 4) pa, Ay = k, ky. 


This result is corroborated by the form of the matrix [G (A | B)]?. G(A | B), 


which is 
Vy F ATA 
Kes | Fe 


1 


ki kz Ey » 
1 2 


v, . BTB 





Corollary 3.2.3—If A = Bin Corollary 3.2.2, [G (A | A)]? represents a SRGDD 
with parameters 


V = 2b, B= Vv, R= vk, K = 2r: A, = ve, A, = ke 


ives Vrtiag ee e 


If Am,, and Bpyq ate any two matric 
G (A/B) is defined if and only ifm + p= 
is defined if and only if 


es over F, then the product G (A | B) 
n + q and the product G(A/B) G (A | B) 


mp = ng. 
If A and B are the incidence matrices 


i of BIBDs, then a necessary and sufficient 
condition for the existence of any or both of the above product— matrices is that 4 and 
B are SBIBDSs. 
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In what follows, unless stated otherwise, A (v,, ki, A,) and B (¥,, Ks, Ao) will be 
assumed to be SBIBDs, 


It has been shown’ that the Kronecker sum of incomplete block designs gives 
ternary designs. We obtain below, an n-ary design as the Kronecker sum of the 
squares of two SBIBDs. 


Theorem 3.3—The product 
G (A | B) G(A/B) 


is the Kronecker sum of A® and B* and is the incidence matrix of a 3-PBnD with 
rectangular association scheme, where n = ki + k, + 1. The parameters of the 
design are 


V = Vy Vo = a R a Vv) ke a Vo k* = K 
A, rs ky ((ky —s 1) Ai + k,) + Vy A, (k; 4+ Ks ap A») + 2k} kS 
Me = v2 Ay (ky + ki — As) + yy (hp — 1) Ay + Ky) + 2ki Ke 


A; = Vo A; (ki ~ ky = Ai) + Vy Ag (ke + ky = As) + 2k} k* 


ny; = Vy — ie nN, = Vv, — l, n3 = (, — 1) (¥, — ty 


In the above theorem, as well as in the ones that follow, we have extended the 
definition of n-ary designs given by Das and Rao to include those designs for which 
incidence matrices have any of n distinct non-negative integers as its elements. The 
entries in the matrix G (A | B) G (A/B) above belong to the set {0, 1, ..., ki + Ko}. 


Corollary 3.3.\—In the above theorem, if A = B, then a 2-PBnD can always 
be constructed, where n = 2k + 1. Its parameters are 


: V=v?= B, R= 2k = K 
Ay = vA (kK? + k — A) + vk* + 2k! 
A, = 2vA(k? + k — A) + 2k! 
n, = 2(v — 1), np = (v — 1). 


The above design has L2-association scheme. 


Agrawal’ has defined binary BBDs with constant block sizes. In the following 
theorem, we obtain an n-ary BBD with unequal block sizes. 


Theorem 3.4—The matrix G (A/B) G(A | B) 


420 G. A. PATWARDHAN AND SHAILAJA SHARMA 
v, A? 

k, k, BE; v 

21 


is the incidence matrix of an n-ary BBD with two different numbers of replications and 
two different block sizes. The parameters of this design are 


ik Be | 
12 

vy, B2 | 

2 











Vo=yv+vy=B 
Vo k, (k, + k,) 
Vy k, (k, + k,) — Ke 


R, K, 


I 
i 


R, 


I 


A, = v3 (ki + ki — A:) + Ke: kK} vo 


I 


A, = A, v? (ki + ky — A.) + kt ke vy 
1 


Avs = ky k (V2 kt oe Vy ke ). 


Here 
n=k, +k, +2. 


Corollary 3.4.1\—In the above theorem, replacing B by A we get a group 
divisible 2-PBnD, where n = k + 2, with the following parameters : 


V=2v= B, R=2vk? = K 
A, = Av? (k? + k — d) 4+ kty, A, = 2vk4 
nm, =v— 1, nm, =», 


Theorem 3.5—If A, B, Cand D are four linked BIBDs with parameters respecti- 


vely 
(%, by, ri, ki, Aa, Bet ST) 3, 4 
such that 
y= b, 
and 
v, = b, 
then 


E=[6(A|B) | [@(CjD)r] 


is a rectangular 3-PBIBD with four unequal block sizes. Its parameters are 


P=" B= 2+), R=, + rb hy + 
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K, = 2k, Ky =v, ky, Kz = vers, Ky = V1 rg 
A, =r, +A, ths + py, As = Ay +r, + psy + ky, 
A, + As + py + py 


% 
l 


fee Meee Vi 1 iy (V5 — 1):(v, 1), 


Corollary 3.5.1—Under the condition that A — B= C = D, the matrix F = ET 
(where E is as given in Theorem 3.5) is a 3-PBIBD with Extended Group divisible 
association scheme and parameters given by 


V = 4, B=vy?, R= ky, K = 4k 
Ai = vA, A, = k*, A, = vk 
nyg= 2(v — 1); n, = 2y, ny = 1, 


Lemma 3.1--Suppose that A, B,C and D are BIBDs with respective sets of 
parameters 


eee head) f) se bo 2. 3.4, 
(i) If N= G(A |B), then NNT = AAT © BBT 
(ii) If N = [G(A/B)]’, then NNT = ATA @ BTB 
(iii) If N = [G (A | B) | [G (A4/B)]*] 

then 
NN? = (AAT + CTC) © (BBT + DTD 

when 
vy, = b, and ve = Jy. 

where © denote the R ronecker sum of matrices. 


The above lemma may be useful in determining C,(NN7) and hence finding 
necessary conditions for the existence of the designs denoted by N. 


4. SUMMARY 


The following table summarizes the results obtained in the last section. 





Sr. 
No. Source Operation Result 
re Thin..3,! G(A | B) 3-PBIBD with rect. asso. Sch. 
and 2 unequal block sizes. 
pd Cor a.).! G (A | A) 2-PBIBD with L, asso. sch. 


(Table continued on p. 422) 
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Sr. No. Source Operation Result 
cf Gores<h2 Submatrix of G (4 | B) Singular GDD 
4 Thm. 3.2 G (A/B) BBD with two replicates. 
a Goresicil G (A/A) SRGDD 
6. Thm, 3.3 G (A | B) G(A/B) = A* @ B? 3-PBnD with rect. asso. sch. 
ico Corad,3:) G (A | A) G(A/A) 2-PBnD with Z.-asso. sch. 
8. Thm. 3.4 G (A/B) G(A | B) n-ary BBD with two replicates 
9. Cor. 3.4.1 G (A/A) G (A | A) GD-PBnD. 
10. Thm. 3.5 {G (A | B) | [G (C/D)]"} 3-PBIBD with rect. asso. sch. 
and 4 block sizes. 
oF Cor. 3.5.1 {G (A | A) | [G (A/A)]™} 3-PBIBD with extended GD 
asso. Sch. 
a a a 
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In this paper we have studied the properties of Finsler space F” of dimension 
n 2 4in which Cij, is of the form Ci. = Gijy) (his an + Ci C; by}. Since 
Cijx of any three dimensional Finsler space is of the above form, any n-dimen- 
sional Finsler space F” (n ? 4) will be called a C3-like Finsler space if its 
Cisx is of this from. 


1. INTRODUCTION 


There are three kinds of torsion tensor in Cartan’s theory of Finsler space F". 
Two of them are (A) /v-torsion tensor C;,« and (vy) hy-torsion tensor Pi; which are 
symmetric in all its indices. It is obvious that F” is Riemannian if the tensor C; jx 
vanishes and Deicke’s theorem shows that F” is also Riemannian even if the tensor 
C; (=g* Ci;x), contracted by the reciprocal tensor g/* of the fundamental tensor gjx, 
vanishes*”"*"!7_ As is well known, in the 2-dimensional Finsler space F*, C;;, is always 
written in the form 


LCi jx = [Mm jm is bbvk) 


where L(x,y) is the fundamental function, m; is unit vector orthogonal to the element 
of support y‘ and / is called the principal scalar by Berwald’’*’’. A Moor studied the 
property of n-dimensional Finsler space F” (n > 3) in which Cy, is of the form 
(1.1)'?""8, In this case, however the third curvature tensor Sjjx vanishes, so that ” is 
reduced to a Riemannian space in virture of Brickell’s theorem provided that the 
metric is symmetric* *, 


In a 3-dimensional Finsler space F*, C;,. is always written in the form 
LG C1 jk = Hmmm — J Giijk) {mm ;nk} 


+ I Gojr {minjng} + J ninyne «»o(1.2) 
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where Gcayjk){ } denote the cyclic permutation of indices /, /, k and addition. H, I 


and J are main scalars and (Jj, my, 1) is Moor’s frame”*”®, Here /) =0; E is the unit 
vector along the element of support, m is the unit vector along Cy i.e. m = GJIC 
where C? = g‘/ C,C, and m is a unit vector orthogonal to the vectors /; and m;. Since 
the angular metric tensor /,, in F* can be written as 


hy; = mm, + nn, atl) 
we may write (1.2) as 
Cijk = Cayk {hij de + CiC iby} pal D4) 


where 


] J 
a, = z (im + “3 (Mk ) 


and 





l H 4J 
bk = 5G8 i( 5 — 1) m +m}. 


Many workers®’!°"!6 have obtained various interesting special forms of C;;x. They 
are 


Generalized C-reducible : 


Cijyk = Cain {Pry Ox. mati 


Quasi C-reducible : 





Cijk = Cayjry {Piy Cx} .--(1.6) 
Semi C-reducible : 
Ci = Cask) i cx + if CC jC, Ba), 
n+] ts 
C- reducible : 
Cijk = a Cayry (hi; Cr} ETB) 
C2-like : 


] 
C1 jk a Cz CC jC; 


where P;, in (1.5) and (1.6) is a symmetric tensor, 
that Pig = Oo = O and P, 


...(1.9) 


‘ Q; in (1.5) is a covariant vector such 
q in (1 7) are scalars satisfying p + q = 1, 
The purpose of this Paper is to study the properties of F” 


of dimension n > 4 in 
Cijk is of the form (1.4). Since Cijx of any three dime 


nsional Finsler space is 


which 
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of the form (1.4), an n-dimensional Finsler space F" (n > 4) will be called a C3-like 
Finsler space if its C;;, is of the form (1.4). 


2. THE PROPERTIES OF C3-LIKE FINSLER SPACES 


Definition 2.1— A Finsler space F" of dimension n > 4 is called a C3-like Finsler 
space if there exist covariant vector fields a, and 5, in F" such that its (h) hy-torsion 
tensor C;;, can be written as 


Cijk = Cyr (hij ae + CiCjbx}. hack) 
Since C;,, is an indicatory tensor it follows that a, and b, are indicatory tensor i.e. 
a, = b) = O. If by isa null vector then contracting (2.1) with g’* and putting b, = 0 
we get a, = ae. Hence it reduces to a C-reducible Finsler space. Thus a C3- 


like Finsler space is a generalization of a C-reducible Finsler space. Again if a, isa 
null vector then contraction of (2.1) with g* gives 





C; = C2 b + 2b; Cc} C; asf2i2) 
; , Ee ; ] : C; 
which after contraction with C‘ gives b, C’ = EE Thus (2.2) gives 5; = 3Cr° 
Hence for a, = 0, the form (2.1) of Ci; reduce to 
| 
Ci = cE CiC Cx. 


Therefore the C3-like Finsler space is also a generalization of C2-like Finsler space. 
Furthermore if we put a, = PC,/n+1, and b, = gC,/3C? then the form (2.1) reduces 
to the form (1.7) i.e. F" becomes a semi C-reducible Finsler space. Thus a C3-like 
Finsler space is also a generalization of semi-C-reducible Finsler space. Contracting 
(2.1) with g** we get 


KC; =(n+1l)a+ Ch ne ge 
C2 

where K=1—b,c’. If bjc’!=1, then K = 0 and hence from (2.3) we get a, = — TEL br. 
Substituting this value of a in (2.1) we get 

Cijk = Cage) {Pry dx} .oo(2.4) 
where 

2 
Pi; = CC, ng en: hij. 


The form (2.4) is the form (1.5). A Finsler space whose Cij.x is of the form (1.5) is 
called generalized C-reducible Finsler space'®. The examples of non-Riemannian 
generalized C-reducible Finsler space is given by Okada and Numata". If 


a = {ar, (x) y'y!Ph 
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and 
B = {bi, (x) yiy pr 


be two independent functions then an n (+2) dimensional Finsler space F" with a 
metric L = L (a, 8) is generalized C-reducile. 


Hence we have the following : 


Theorem 2.1—Every C3-like Finsler space is a generalized C-reducible Finsler 
space provided b; Ct = 1. 


If b; Ct # 1 then from (2.1) and (2.3) we get 








Cijk = Cajry {Qty ax} + diddy asi as 
where 
l 
Oy tee ig 8 (tS ir (a; bj + a, bi) 
2(n+ 1)C? 
+ BAO nb, (2.6) 
and 





4 \1]3 


. Now we shall give below an example of a Finsler space whose (h) hy-torsion tensor 
Cijx is of the form (2.5). Let F" be a C-reducible Finsler space with metric function 
L (x, y). If F*" be another Finsler space whose metric function L* (x, y) is given 
by 


L*9 (x, y) =L3 (x, y) + (b; y')3 


where 5; (x) is one form in F". The (A) hy-torsion tensor rales of F*" and those of F" 
are related by 


x Pq? . 
Cin = P Ci + SE (hs me + hig my + hyx my) 
2p3 — | 
FP mim ym, al YF 3 
where 
oe 7 ye 

m= qi, Phu,gm 7, q— 5 pay 

Thus putting 


l 
Cie ye : 
tah (GN BET, (11; Ce + ix Cy + Ay C)) in (2.7) 
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we get 
Cae = Case (Qi; a} + did jde 
where 
O1, = phij, ay = a + Lm and d; = — (Pe Pp Jim, 


A Randers space is a Finsler space whose metric funciion L (x, y) is of the 
form 


L (x, Y) = (ar; (x) iy)? + fi (2) yf 
where a;,; (x) is a Riemannian metric tensor and fj (x) is one form. 


A Kropina space is a Finsler space in which 


L(y) = Udy 


It is well known*"’* that the Randers space and the Kropina space are C-reducible 
Finsler spaces. Thus if the metric function L* (x, y) of F*" is given by 


E** (x, y) = {(arj (x) yy? + fix) YP + (Cr vf 
or 


*3 Pit (aij ix yyy wrd\S 
ee Choe. 62) 


then the (A) /y-torsion tensor of F*" is of the form (2.5). 
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Some approximate solutions of the Duffing equation is considered. The 
approximate solutions are exact for a particular initial value. 


1. INTRODUCTION 


Many physical phenomena are modelled by ordinary nonlinear differential 
equations. In order to analyse the behaviour of the physical systems, we need the 
solutions of the differential equations which is extremely difficult in the nonlinear 
case. However there are methods to obtain approximate solutions of nonlinear 
equations. Peters‘ and Usher* have obtained some approximate solutions to the 
equation 


xturx= —ax? Al) 
subject to the initial condition 
x (0) = A, z(0) = 0 ...(2) 


But their approximate solution is not equal to 
x(t) = —o*/a when A = —*/a, 


Such a difficulty has been eliminated by Shidfar and Sadeghi? by using a different 
method. 


In this paper we shall obtain some approximate solutions of the Duffing 
equation 


x + wx = BPxt ...(3) 
with initial conditions x (0) = A and « (0) = 0 
by suitably adopting the method discussed in Shidfer and Sadeghi.” 


In pursuit of this we put 


x =u —o/V 38. (4) 
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Inserting (4) into (2) and (3) we obtain 
it JV 308 u2 — B72 u3 = 2 w3/3 V/ 38 


u(o)= A nea Wolo eo 


Next we shall discuss the series solution of equation (5). 


2. Series SOLUTION 
To solve equation (5) let us investigate the solution in the form 
u(t) = co +c) Sin wt + c, sin? wt +... 
where the coefficients Cy, c,, Cc, ... are constants and to be determined. 
Substituting (7) into u (0) = A + w/V 38 we get 
Co = A+ fv 38. 


to} 
.-.(6) 


mats, 


...(8) 


Using (7) into LHS of (5) and equating the sum of the coefficients of each sin” wt, 


n = 0, |, 2, ... to zero, we have 


20? cz + V 308 co? — Be? = 2 w/34/38, n= 0 


i et a hel 
Cn+2 (n + 1) (n = 2) Ch a(n + 1)(n + 2) b, 
B2 


= wat lint) (Doc, +... + b,co)n > 1 


where b, = Coc, + ... + CpCo by = cé. 


From u (0) = 0 and (7) we can immediately concluded that c,; = 0. 


The second equality of (9) and c, = Oyields cs = 0,c,; = 0... successively. 


the coefficients with even subscript we obtain 
Co = A + o//38 


roe j 
OA ae ed 


ey ea AB 2 2/a2 
©= Foe (4? — w/B2) (A? + 0/82), ete, 
This gives the solution 


A B* 


xX()= At ee (4? — w?/B*) sin? wt + 





for equation (3) subject to the initial condition (2), 


...(9) 


For 


...(10) 


(1) 
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Remark 1: As we see the coefficients c,, cy, vanish for A = + /B and it can 
be readily shown that the other coefficients also vanish for this value of A. 


Therefore x (t) = + /8 for A = + /8. Hence the solution is exact for this 
particular initial conditions. 


3. CONVERGENCE OF THE SOLUTION 
It is interesting to note that the series solution (7) convergent for allt. Infact, 
in the case cp >0, c, >0 and 8 = — y, we obtain from (9)c, >0forn = 4,6,... . 
By the second equality of (9) we can immediately conclude that 











0co co oo 
S ervernen= F wa rt(S o) 
n=2 nw2 n=2 
2 ~ 3 3 2 
ge) coo a a 
a=0 
(12) 
or 
“ 2 2 =< 3 
mt (S JS s) 
n=0 n=0 n=( 
2 
= Sasa Co" + a Pt —s Ae 


co) ...(13) 
Therefore, in this case, the series of coefficients of series solution (9) converges. Now 


if we put 


1 
Cr mile, c = 0, C=C, 


/3 8b} 
@ (n + 1) (n +2) 


1 n? 


1 
n+2 (n + 1) (n + 2) Fag 


aid ne bt Cc 


B? 1 
+ GF wea % & oy a SI 


on 


where 


n 0’ 


then the series = c. converges. Since |c, | < c} it follows that the series solution 
n=0 

(7) is absolutely convergent for all t and consequently the series solution (7) converges 

for all t. 
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4. Power SERIES SOLUTION 


Finally we seek a power series solution in the form of power series say 


u(t) = = 7 ae fe eet 15) 


n=0 
for equation (5) subject to the initial condition (6). 


Inserting (15) into (5) and using the initial conditions (6) we can conclude 


a;.=.0, ag 0a, 0 2 32(16) 
and 
a = A+ w/V3 8B 
2 a, + V3 wB as — Bah = 2 w5/3 +/3 B. 
= 9 2 
ae ais SE a CESSES da, Ara 
+da)n=2k,k>1 Ate) 
where 


dy, = Apdn + 2 Qn-g +... + Gn Qo. 


The coefficients (17) are vanished for A = -+ «/8 and the correspending solution of 
eqn. (3) becomes 


x(t) = + off. 


Remark 2: The solutions (7) and (15) are identical. This can be easily verified 


by substituting the Maclaurent series of sin wf into (11) and writing it asa power 
series: 
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GENERALIZED ORTHOGONALITY RELATION FOR THE 
FLEXURE OF SECTORIAL PLATES 


H. SRINIVASA RAo 


Department of Mathematics, Karnatka Regional Engineering College 
Srinivasnagar 574157 


(Received 2 January 1987; after revision 29 June 1987) 


In this paper using some refined plate theories a generalized orthogonality 
relation for the flexure of a sectorial plate is derived. 


INTRODUCTION 


Eigenfunctions of the biharmonic equation and orthogonality or biorthogonality 
relations for satisfying the boundary conditions have been extensively used in the two- 
dimensional elasticity theory and in the flexure of thin elastic plates. Smith’ has pre- 
sented an orthogonality relation for the biharmonic eigenfunctions associated with a 
rectangular strip whose long edges are stress free. Greenberg? has brought out some 
orthogonality relations satisfied by the biharmonic eigenfunctions associated with a rec- 
tangle whose long edges are stress free and a sector plate whose radial edges are stress 
free. Prokopov® has shown that the orthogonality relations obtained by Greenberg exist 
even ifthe long edges of the rectangle are built-in. A very general method of develop- 
ing a biorthogonality relation for obtaining the solutions for a semi-infinite strip whose 
long edges are stress free is provided by Johnson and Little*. Little and Childs’, Klemn 
and Little’ and Klemn and Fernandes’ have successfully used the method of biorthogo- 
nality relations to provide solutions of various boundary value problems for different 
geometries. And newer orthogonality relations for the boundary value problems 
associated with a rectangle or a cylinder or a curved beam have been discovered by 
Nuller®, Kostarev and Prokopov’, Fama’? and a few others. 


A study to develop eigenfunction solutions has not been made for the plate 
theory which takes into account shear deformation of the plate. Further, the available 
literature does not show the existence Of generalized orthogonality relations satisfied 
by these functions. These relations are indeed very elegant and can be directly used 
to reduce boundary value problems to a system of linear algebraic equations. Accor- 
dingly the theme of eigenvalues, eigenfunctions and eigenfunction expansions in relation 
to the refined theory of plates has been explored in considerable detail in the doctoral 
thesis of Rao". In what follows, a derivation of the generalized orthogonality relation 
for the flexure of a flat curved plate taking shear deformation into account 1s presented. 
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Basic RESULTS 


The basic governing differential equation in polar coordinates for the flexure of 
plates in the refined theories of Reissner!*, Hencky’* and Mindlin™ are 





























aw = £2 (i =" ) 5 ada ...(1a) 
Ay = a. d ...(1b) 
where 
= — ge ote PDO Q, vasteel 
PPE ee ~ a aDaSD ap ...(2b) 
mono [Bee (LE ES) 
Got (CF Sp) ee ee 
mr Le oot — Be + 200) 09 
= - DF (m+ t Se Milan i ...(2f) 
Qo=—D 7 say - Eo dan og ...(2g) 


Here »? is the shear rigidity factor introduced by Mindlin. For Reissner’s theory we 
a m=2,n = Hi2 and 7? = 5/6. For Hencky’s theory we set m = 1,n =O and 
n° = 1. For Mindlin’s theory we set m = 1, nm = O and ny = 5/6. 


FORMULATION OF THE PROBLEM 

We now assume that the curved 
0<6<58 is free from the load q (r 
(1) and (2). Defining a new function j 


plate occupying the region a <r<b, 
6) and accordingly set g = 0 in the equations 


P= Aw 
a3) 
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the plate equation (la) may be written as 


Ap = 0 ...(4) 


(1b), (3) and (4) provide three partial differential equations of the second order. Alter- 
natively we want to formulate three partial differential equations involving either the 
variables w, Mz and ¢, or the variables Qo, «2 and M,s. To this end we connect p and 
with the above mentioned variables. It is shown in Appendix (A.1 — A.2) that 


M 
p= se — (1 = 'p) € et py 


and 


v a ek Mre 


D on Se ...(5b) 


where the prime (’) denotes differentiation with respect to r. From (2a) 


ow h? 
A(retr-)= arpaam 80 


Using the operator relation 
0 7) 
a(r-s-)=(" =) a +2 


) and the relation (2f) to replace A (r Q,) we get 





é 
to replace the term A ( i 


or 
A Raed Se a ur ( phase 
A (r ¢,) ina ane “| 2p he Fae Bil aes )=0. ...(6) 
From (2g) 
op ou 
7 ON hove car -emallaml ey a 


ra 
and remembering that A (=) = 0 we have 
12%? av ) a 7 
A(r Qe) + 2 ( r5 + 29) =0. A(T) 
From (2b) 
ow h? 
aay af — St aap o | 


op sia aN r 
~ ae ra aa le a 
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Using (2g) and (7) we get 
r ow h? 
sred=5(O+ FF )— arpa =p) 


12? (ay ) 
h? ( ér bee, 





ur 9 I+pe r a a 
A (ra) =r + Le D "Or Tn D = 0, ...(8) 


Thus the differential equations (3), (4) and (6) are in the variables w, «and Ms, The 
equations (7), (8) and (1b) are in the variables Qo,€o and Mo. 


We now seek the eigenfunction expansions of the functions w, p, b, €r, €9, M,,Me, 
Mo, Q, and Q¢ in the form 


| w, P, €ry €6 | = am exp (Ax 6) | We Pr, Erk, coe | -(9) 


and 


2 ral v, Mr, Mer Mro, Or, Ov | 


a 7 ak exp (Ax 6) E Mk, Mok, Mr ox, On, On | ...(10) 


where A, is an eigenvalue. 


Using the above expansions in (Sa), (Sb), (3), (4), (6), (7), (8) and (1b) we 
obtain 


Pye [ Ma + (1 — #) ce | .. (11) 
= — [Mt (1p) ed (12) 
Pew. +trwe + Ay We = r? py ...(13a) 
r Pe trp + A py = 0 ...(13b) 








2 e ‘ 127? r2 1 
feta Oi Fi ] — i )ex— jt r? pi 


12 . 
+ 2rp, — re row + Ay fn =0 ---(13¢) 
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127? 
ht 


( vi + 2r de ) sat) ...(14a) 


r? OF, + 3rQi, + Qu +A, Oo + 


2 
Se gta Car thee AS tog — 6? ox 


Lr 


ut ; 
+ eee + 





4 
TR loge ...(14b) 





th oe ok tak, Pt a d, be = a r? Wk. ...(14¢) 
To pose (13) and (14) as system of adjoint differential equations we introduce 
Ue = Wr, Ux = Pr = =| Ma +(l1—p €, ] » Use = érx 
and 
Vik = Qon, Vox = €ox, Vox = Ve = — | Mrox + (1 -- ») ‘esl F 
(13) may then be written as 


r? UN, +r U;, —ru,= — At Usk vs-(15a) 
rig + rU,, = — Ai Ux ...(15b) 


r?U*, + 3r Us, + ( 1 — 12y? 4, =~) Use 


] te 127? 











at ree Oy = At OF, © ...(150) 
(14) may be written as 
aye PT pe Ve th Or Voy | — A Ven --.C 162) 
f V5 + 3r Vie V. ie iF Tigi rvs} = k ' 1k via 
+ Ree ers 

Beer et Ve = OV po 

4+ . 

+5 or ea s M Vox ...(16b) 
2 

pe be Vig vay Vaz = — Xe Vox. ...(16c) 
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THE ORTHOGONALITY RELATION 


Equations (15) and (16) may be written in the vector form 





L(O) +02 T =0 ia) 
and 
L*(V,) + BM. = 0 (17b) 
where 
{ ) f } 
| Renee a U,-r Ur | 
Uk | 2 “ | 
. | fall! SAE ead #5 
| |? 127? 
ve | 2 US, + 3n.U; .+(1 Sipe ) Us. 
| | | 
! | | I+p ne aes 
i l +2r U4. i i—p, 3 Cpa hi —_—r? G., J 
and 
f } ( - 
| | Me Spe er = (Hv; : +2rVax) } 
| V | | 
| Va Wee : | 
¥ : | +3 Vy, + Kea — 8? Van 
ha A | L*(V,) = | ; | 
. +p Sa 
| Va | i+(j Tau) (1 Vin + 2r Pan) + 2r Pa 
| | | | 
| : ; 12y* 2 
[ j Rane: a +r ke. — Ie rf V3x j 
It follows from (17) that | 
b 
6 


[tz (Ux). V; — L*(V,). 0] dr = ( Ai — 23 | (Ox.V ;) dr. 
Expanding the integrand in full and carefully integrating by parts (18) reduces to 


b 
( rj ad Xi ) J (Wk Qo, — €oy Mox — €rx M,0,;) dr 


=| p32 = ‘ 
[+ (Uj, Vij Vix ES + Un Voy — Uy Vis ant roe Vs, — Ux 





; ag 1 
+ r (Us, V3, U2, Vi; — Uy Vi;) cet = : (r? Un, V3) 


...(18) 


eet 


(equation continued on p. 439) 
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12 2 b ° 9 
7 (r? Uix V3;) [+( Aj — Ay Ja — p) EE Eo; iz .. (19) 


The components of the vectors U,V, U’, V’ are related to the functions We, €rk, €ox. 
Mok, Mrox, Qox and Q,x as detailed in the Appendix (A.14 — A.25). 





Upon considerable simplification of the right-hand side of (19) making use of the 
above relations we obtain ultimately 


b 
( 2 re dy | (We Oey — €6; Mox — €x Mraj) dr 


= r A; (Wy Ore — We Qrj) + 1 Ak (€0; Mrox — €0% Mroj) 
b 
— r Aj (€r; M,. = Cphp M,)) | . ...(20) 


The following combinations of homogeneous boundary conditions may be prescribed 
on the curved boundaries of the plate ; 

Case (i): Clamped: w = «, = «, = 0 

Case (ii) : Simply supported : w = M, = Mro = 0 

Case (iii) : Free: OQ, = Mr = Mro = 0 


Case (iv) : Elastically clamped : w = ¢9 = 0, M, = we, where « is the moment re- 
quired to produce a unit rotation «,. 


Case (v) : Elastically supported: M, = Mre = 0, QO, = Bw where 8 is the shear 
force Q, required to produce a unit deflection w. 


The right-hand side of (20) vanishes for all these conditions giving the desired ortho- 
gonality relation 


6 
| (We Qo; seg SF, Mex — €rk Mro)) dr = 0 for all 7 = k. .(21) 


Interchanging j and k above and adding the resulting equation to (21) we get the gene- 
ralized orthogonality relation 


b 
| RO ares hin: Me, 
a 


+ wy Oox — €0x Mo; — rj Mrox) dr = 0 .+2(22) 


valid for all j # k. 
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The application of these orthogonality relations to boundary value problems in plate 
theory will form the basis of subsequent work by the author. 


REFERENCES 


R. C. T. Smith, Aust. J. Sci. Res. 5 (1952), 225. 

G. A. Greenberg, J. appl. Math. Mech. 17 (1953), 211. 

V. K. Prokopov, J. appl. Math. Mech. 28 (1964), 351. 

M. W. Johnson, and R. W. Little, Quart. appl. Math. 22 (1965), 335. 
R. W. Little, and S. B. Childs, Quart. appl. Math., 25 (1967), 262. 


SNAAWR YN 


J. L. Klemm, and R. W. Little, Siam J. appl. Math. 19 (1970), 712. 
J. L. Klemm, and R. Fernandes, J. appl. Mech. 43 (1976), 50. 
B. M. Nuller, J. appl. Math. Mech. 33 (1969), 364. 
9. A. V. Kostarev, and V. K. Prokopov, J. appl. Math. Mech. 34 (1970), 902. 
10. M.E.D. Fama, Quart. J. Mech. appl. Math. 25 (1972), 479, 
11. H. Srinivasa Rao, Ph. D. Thesis, I.I.T., Bombay, 1981. 
12. E. Reissner, J. appl. Mech, 12 (1945), 69. 
13. H. Hencky, Ingr.-Arch., 16 (1947) 72. 
14. R. D. Mindlin, J. appl. Mech. 18 (1957), 31, 


APPENDIX 


In deriving the results below expansions (9) and (10) have been used wherever 
necessary 


Using (2a) and (2d) and remembering that 
7) r r 
2, 5 Or 4 1 20, 








ér r r ON 4 
we get 
Me ae 
p= D ons (1 a pu) ‘. .-.(A,1) 
Similarly using (2f) and (2g) we get 
004 l a see 
or. r- ws, oe : 
7? 
= — 12 Re } using (1b). 
From (2b), (2e) and above relation we get 
cz a Mo ‘ 
D >> — Cee ...(A.2) 
Using (2a), (2b) and (2e) we have 
Mz 








=k 1 — A(t de, €6 ; 
D 2 r 08 a r 7 €o ) ° Aas) 
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From the above relation it follows that 


‘ _—s 2 Myo ae Ce, 
Der eee + ~-(e — 3). 


Eliminating ¢, between (A.2) and (A.4) we get 








D D r 
From (2a), (2b) and (2d) we get 


meee rece 5). 











Ms €, 1 deo Ce, 

D oe r 7 r 06 = or 
Using (2c), (2d) and (3) we get 

M, M 

D - Ke = — (1+ 4)>?. 


Eliminating p between (A.1) and the above relation we get 


Me M, 2 ‘ 





Eliminating M, between (A.6) and the above relation we get 
ee ee Oo 
ange PD ee ane 
From (A.8) and (A.9) we get 


M M, é 
>= Sere Ci =P) («+ 3). 








From (A.7) and the above relation we get 


M, O€g 
p=-  —a—w(at 3 ). 








Using (2f ) and remembering that p = Aw we get 
1 1 ay 
PE Oat ad 


Eliminating Qe between (2g) and (2b) we get 


cia amiepeenitae (Todt. o 
eer ae a (6atXl =p) \ r08 ar 
Lc. 
ot (lw, 1 ow) Lo 
1 lo eee meen | al eT r 00" 


) 
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...(A.4) 


.. (A.5) 


..(A.6) 


Sel (Win 


.(A.8) 


..(A.9) 


..(A.10) 


.(A.11) 


..(A,12) 


. (A.13) 
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Using the definitions of the components of the vectors U and V we get 
U.~ = We 
Usk = Pk = — Mr —(1 — p) (ere + Ab €ox) using (A.11) 
Us. = erk 


2(1 — us A : 
Vir = Qo = be et) €ok + on Wk ) using (2b) 


Vox = € 6% 


| = 
Vok = Ue = Mr + — P (At eck — ox) using (A.5). 





Using (2a) we get © 


he 


i? = Wa.= = Cpl a tthe wr Or. 


10k 


Using (A.12) and (A.5) we get 


; , A, M, — 
ase =w=-O,.-— a Malt | edt 3% 
r r 


(Ax €rk — € 9x). 


From (A.9) we have 


‘ rte ‘ = a 
Uy — ©, = Ma — “p (erk + Ak € 0k). 


Using (A.17), (A.4) and (A.20) we get 





Yh ee ‘ wy 67* (1 — 
Pe Ee Tt See ee ed 


a Bell) See 


r r 
Using (A.4) 
OF abe 2 ] 
VV. = “<< % a Mrox + 5h <ok — Ak tte). 
From (A.13) and (A.11) we get 
oe 6y2 (1 — 
Vra=y= (en 20 ae Wk ) 


+ [ Mn + (1 = B) (ere + Ag eax) i: 


...(A.14) 
AAS) 


..- (A.16) 


...(A.17) 
wo(A.18) 


avst{Aslo) 


...(A.20) 


..(A.21) 


...(A.22) 


..(A.23) 


..(A.24) 


--.(A.25) 
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A NOTE ON THE CHANDRASEKHAR’S X-FUNCTION 


N. SUKAVANAM 


Department of Aerospace Engineering, Indian Institute of Science, Bangalore 560 012 
(Received 12 March 1987) 
A differentiability property of the Chandrasekhar’s X-Function is analyzed. 
INTRODUCTION 


Chandrasekhar’s X and Y Function play a vital role in the theory of radiative 
transfer’ and neutron transport theory”. These functions are defined as follows : 


lim \ 
X()= og J 2) | 
and : r nla 
x 1m | 
VAAL Fs ccbl, 7 ay: Cx. 8) j 


where J (x, z) is the unique solution of the integral equation 
t 
J (x, z) = exp (ixz) + J K (x—y) J(y, z) dy. a Alaa) 
0 
We denote the functions ¥(z), Y(z) and J(x,z) by X(z,t), Y (z, t) and 


J (x, z, t), respectively, to emphasize their dependence on ft. In this note we analyse 
differentiability property of the X-function with respect to f. 


Let 7; be an integral operator from L, [0, ¢] into itself defined as 


t 
Ti (f) = | K(x — y) f(y) dy. 
The adjoint operator T* defined by 


T (f= j K@ —*) £0) dy). 


Assuming that 7; is real and symmetric, i. e. 
K (x) is real and K (x) = K (— x) (1:3) 


K (x) is locally square integrable ...(1.4) 
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K (x) E li (- oo, °°) ge) | 
and 


{' —.K (a) 56.0, cee (1.6) 
where K (@) denotes the Fourier transform of K (x). Vittal Rao*® proved that 
x0) =a, X01 (1.7) 
where « (x, t) is the unique solution of the equation 
a (x, 1) = K(x) + K(x ohn Gea) ee (1.8) 
In this note we generalize this result for a larger class of operators, which 


satisfy, instead of (1.3) and (1.6) the following conditions 


T; is normal for all t ret Leal 


bei) K (@) || G0 00 ities: wl Gy 
The generalization is established by the following theorem. 


Theorem 1—The unique solution « (x, t) of (1.8) is differentiable with respect to 
t for almost every ¢ and its derivative is a (t, t) a: (x = ¢, t) in the following sense. 


Lt 
call a(t, )aq@—4n| dx =0 
0 ...(1.9) 
where «, (x, ft) is the unique solution of the integral equation 


a (x,t + h) — «a (x, #) 
a ee 





11) = K(X) +. K(x — ») a (41) aby. (1.10) 


Theorem 2—For real z 


4] ee | 
“or * ©, t) = exp (itz) a* (t, 1) X* (z, t) dah) 
a SE, 
oF X* (z, t) = exp (—ilz) « (t, ty Xz 2) ...(1.12) 


where X* (z, th= Lt é 
| x eae J* (x,z,t) , J* (x, z, t) is the unique solution of the equa- 
tion 


J* (x, 2, 1) = exp (— ixz) 4 [KO ~x IG, 20 dy (1.13) 
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and «* (x, t) is the adjoint of «* (x, 1) given by 


a* (x,y) = K( —x) + i KO — x) a* (y, t) dy. ...(1.14) 


PROOF OF THE THEOREMS 


First we note that | is not an eigenvalue of the operator 7 and hence (1 — 7;)™ 
exists. 


Proof of the Theorem 1—From (1.10) we get 


a(x —t,1)= K(x —N +) Ke —t—y) m(y,1 ay 


eee ae Kies) 1,0) ds. vee(2.1) 
Hence 
a, (x — t,t) = (1 — T)1 (K(« — 2). (2.2) 
Now 
t_h 
a — « _ a(y,t)—a(y,t—/ 
(x, t) (x, t — h) =(K@-») (y, 0) pO! 1) 
+ = K(x — y) «(y, t) dy 
h : 
i-h 
Hence 
A she = (1 — T-n)7} i | K(x — y)«(y, 1) dy |. 
...(2.3) 
Define 


— a(t, t) a(x — f, t). 


nr ee ; SE ae 


Using (2.2) and (2.3) it can be easily verified as in Theorem 1 of Vittal Rao® that 


t 
a al 
1 AO 


a (x, t) — « (x, t —h) 
h 





h ee | Da) 








— a(t, t) a, (x — t,t) | dx = 0. 
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Similarly (1.9) can also be proved. This completes the proof. 
Proof of the Theorem 2—Let RF be the resolvent of the operator 7*. Then 
‘ 
R* (x,y) - FQ —® = | R* (x,s) Ki —s) ds. (2.4) 
0 


From (2.4) an (1.8) it follows that 


R* (0, y)= a(y,f) . ...(2.5) 


The solution J* (x, z, t) of eqn. (1.13) is given by 


f 
J* (x, z, t) = exp (—ix2)+ | R* (x, ¥) exp (— iyz) dy, ...(2.6) 
0 
Therefore 
X¥*(z,t) = Pee J* (x, z, t) 
t 
ger Sale | “(y,t) exp (—iyz) dy. (From (2.5)) 
0 
Similarly 
ra 
eG 2, t = 1 * , ; 
(z, t) + J a* (y, t) exp (iyz) dy. S(27) 
Now 


fen ge ee 
X* (Z, t) — X* (z,t — h) 7 j a (Ys 1) — a (yt — A) 
h ae 
tt) 


t 


eee i ’ u 
exp (—iyz) dy + r | «(y, t) exp (—iyz) dy. -+.(2.8) 


From T 
m Theorem 1 it follows that for real z the first term j in the right-hand side of (2.8) 
converges to « 


t,t — 
(1, 1) ) j a1 (y — f, t) exp (—iyz) dy and the second term converges to 


% (1, t) exp (—itz) as h — 0. 
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Hence 


‘ 


a X* (z, 1) = a (f, t) [exp (—itz) + | a (y — t, t) exp (—iyz) dy | 


Substituting »y = t — x and observing that 


a(— x, f) = a® (x, 2) 


we get 
a ’ ——_ : ; 
oa hy X* (z, t) = a(t, t) exp (—itz) | 1+ | a* (x, t) exp (ixz) dx ; : 


This gives the formula (1.12). Similarly (1.11) can be proved. 


A relation between the X-function and the eigenvalues had been obtained for 
symmetric integral operators*. Many spectral properties of symmetric operators had 
been extended for normal operators by Vittal Rao and and Sukavanam‘. The results 
of this paper may be of use in obtaining a relation between eigenvalues and the 
X-function for normal operators. 
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BESSEL TYPE FUNCTIONS INVOLVING DIRICHLET SERIES 


A. G, Das AND B. K. LAHIRI 
Department of Mathematics, University of Kalyani, Kalyani 741235, West Bengal 


(Received 23 March 1987; after revision 19 October 1987) 


In this paper we construct an ordinary linear differential equation in the com- 
plex plane whose solutions in terms of Dirichlet series, found by the method 
developed in Das and Lahiri [Rend, del Circolo Mat. di Palermo, 33 (1984), 
425-35], appear to satisfy certain relations analogous to the Bessel functions. 


1. INTRODUCTION 


In an earlier paper® we considered a class of ordinary linear differential equations 
in the complex plane for which we obtained solutions by a process in terms of Dirichlet 
series convergent in a half plane. For the analysis of that process we used arguments 
by Cesari’”’. In the present paper we take into consideration a particular second order 
ordinary linear differential equation (4) in the complex plane for which first we exhibit 
the solution in terms of Dirichlet series which results from the process by Das and 
Lahiri*. The same equation (4) could also be solved by power series and the correspon- 
ding solutions would have certain interesting properties. Furthermore the solutions of 
(4) in terms of Dirichlet series could be converted to power series by suitable transfor- 
mations. But our main interest in the present paper is to formulate, independently of 
the power series, some basic properties of Dirichlet series solutions of (4) and to point 
out that these properties are similar to the classical properties enjoyed by Bessel func- 


tions. The present paper may be considered as a direct application of the theory 
initiated in Das and Lahiri’, 


Consider an ordinary linear differential equation in the complex plane 


nw dp 
Tn Ct Ps (2) Gear + + Pa(z)w=0 aeol F) 


where p; (z) possesses, for each i, an absolutely convergent Dirichlet series (ordinary) 
expansion 


co 
pi(z) = &, Puslke 


for Riz > 8 (say) andi = Ae Saray 


The transformation 
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= _ dw d"—w 
Dt Cees a) 4 (ieee dz” 9 Va = dz" 


reduces (1) to the system 


d 
= = A(z)» RAG) 


where y = ();, Vo, -.-» Ya)? and 


fo | 0 i 01 
1 Oo Oo 2 0 | x 
4a—=t eg) oles Alc ay, 
| 0 0 0 “ 1 | 
en apy, | 


In an earlier paper® we solved the system (2) in terms of Dirichlet series with a 
general n x n matrix A (z) where we supposed that the entries aj; (z), 1,7 = 1, 2, ..., 


n, of the matrix A(z) possess Dirichlet series expansions each absolutely convergent for 
Riz > 6. 


If S denotes the set of all vector functions w (z) = (Wi, Wo, -. , W,)7 


where 


ioe) 
w(z)= ZX Crk? 
k=l 
is absolutely convergent for R/z > 54, 
co 
w(z) = 2 CaleZ, Ch = (Ck, Cus, +++) Chn)™ 
=1 
then it has been shown in Das and Lahiri® that S is a Banach space with 
co oo n 
wl = 2 feel fed = (2D | ex; | ed. 
k=l kul jul 
The following results are obtained in Das and Lahiri’. 


co 
Lemma A—For any scalar A, there is w (z) = z cx{xz in S such that y = w e# 


satisfies the system 


(I oo — Aly = (Ml = A) 1 eM 


where / is the m x n identity matrix. 
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Theorem A—For every eigenvalue A of A; the system (2) possesses a solution. 


te ae 
The process involves the recursive relations for the determination of the c;’s as 
follows : 


(Al — A,) c, = 0; 


[Ai — (A + log 4) I]e. = — A, ¢;; 

[4, — (A+ log 4)/]cs = — As ¢;3 

[41 — (A + log }) I]ca = — Aye, —A. 63; 

[Ar — (A + log —) Te, = — >) Ames .-.(3) 


where the summation on the right of the nth relation is extended over all pe Bye To 
that r,, s, are factors of n and r,. s, = Ny Atk e=1), Pag Se Rom 2. sak 


In each of the cases when the eigen values 
(i) do not differ by zero or by a logarithm of a positive integer; 
(ii) do not differ by a logarithm of an integer (> 1); 
(iii) differ by a logarithm of a positive integer (> 1), 
we obtain in Das and Lahiri? two independent solutions of (1). 
It was observed in Das and Lahiri? that if 
Y= (Vay Vay vey Vn)” = OM (Wy, Wey orry We) 


is a solution of (2) then y, = ez W, (z) is a solution of (1). 


2. DIFFERENTIAL EQUATION AND SOLUTIONS 


Consider the following ordinary linear differential equation in the complex plane 


d’w dw 
dz? 4 Prue {a? ez — } (n?a? = 1)} w =0 ...(4) 
where « = log } and n is an arbitrary constant, 
The transformation wu—=wW,y, = reduces (4) to the system 


dy 
ae ede) y ...(5) 
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where y = 1 y>)? and 


fapuee Leeda - }) iN ie ibe 


Ai. LA Ay. 2%, gays 


Eigen values of A, are obtained from det (A; — A/) = 0, i. e. from the equation 
AAFD) —(F =H (F+H)-0 


which has the solutions 


If y Fa (11, ye)T 


= er: (v, (z), We (z))7 


oo co 
= e§F (3% crrlez, ZS Cho 4z)T 
k=1 k= 


fore) 
=e > Ch/| KZ 
k=) 


where cx = (Cx1, Ck2)”, is a formal solution of (5) then we obtain using (3), the follow- 


ing relations 

(AI — Ai) c, = 0; 

[4, — (A + log 3) I] c, = — A, 1; 

[A, — (A + log ir)J] car = — A, cor — 1, r = 2, 3, eed, 
since A, = 0 fork > 2. 


For any eigenvalue A, we have 
— (log $)* car--1, 








Cory =—_ = 
(A+ logt, — 5 log} + 4) + log tr + 5 log t + 2) 
a8 ae 
(—1)’ (log 3)" ec; 
: n 
(\ + plog} — log $+4)(A +p log $ + = log} + 4) 
p=i 
If A, = > += ~ log + — 4 and 7 is not an integer, then the actual 
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computation yields 


—Cy me (- ly Cy 
oF fal)! <8 aaa ee ee 
and in general 
(~ I en pe 12 
Cr = ar re een Wr RE moray pa Fe how 2 be sul 
¥ 1.2...r(n + 1I)(a+2)...n+17r) ’- 
; 1 
From the first relation of (6), we have c: = ¢,, (1,A)7 and so taking c; = Prantl” 
we obtain 
oo ( 1) 
n iz —z)r 
nmen (Fosi- i: S resyreaee @ 
r=0 
= e-2/2 Q-nzh See (2-7) 
~, Pr+)Prat+rt+h 
r=(0 
n 
as a particular solution of (4) corresponding to the eigen value A, = 5 log 3} — 3. 
n 
A second solution of (4) Corresponding to the eigen value A, = — — log 4 — 3 


2 
where n is not an integer, is given by 


00 

= e77i2 Jnz/2 (alk eer —z\r 

panimP Sacre cer re CY 
r=0 


Ifn = + m, where mis a Positive integer, then the two eigen values of A, are 
A, = J log 4 — 4 and A, =- S log 4 = $680 that A, +A; = log 3,. 


A particular solution of (4) corresponding to A, = > log 4 — 4 is given by 


co 


yy = e72/2 D-mz/2 oest 98 -z)r 
DS hee > 77 TeaEr ie = 


To obtain a solution of (4) corresponding to the eigen value Az 
— 4, we consider the formal series 


m 
panes. log 4 
co 
JH{A—A)¢+ > s Cklkz} eAz 
k—=2 


satisfying the relation (cf. Lemma A) 


d 
Caz A= Os 4) a Aas 
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The recursive relations in (6) yield co, c 9 » «++» Com—-) With A — A, as a factor. We 
2 


. . * .- 1 
recall that c, = 0 if r 4 2* for some positive integer s. Thus for \ = A, we have c., 
gto eo 


«+ = Com_, = 0. In fact, forI1 <r<m, 

1 1 

ee (—1)" (log 4)" (A — A.) ey, 
(I (A+p log 4 — m/2 log $ + 4) (A+p log 4} +m/2 log 1/2+ 4} 
p= 


Since A + mlog 4 — * log $+ 4 is the only factor in the denominator of cm, 
which vanishes for A = A.,, it follows that 


com, = (=; 13" (log 4)" 61 








TLQ+p log $—m/2 log $+4) Ml A+p log $+ m/2 log +3) 
= p= 


Te) 
is a rational function of A in the neighbourhood of A, and A, is not a pole of cy’. 
Hence c, m + r;,r > 0, can be calculated for \ = A». In general, for r > 1 


—|])mtr 1)2(m+r) 
cm+r, = (=I) (log 3) Ci) 
Tl (A+p log }—m/2 log }+4) TI (A+plogt-+ m/2 log $+4) 
po=lsp#m p™1 


> (~— 1) (log3)*" com, 


I {A-+(m-+p)log} —m|2log} + 3} {A+(m+ p) log} +-m/2 logt +4} 
p=) 


...(10) 
Choosing ¢;,; = (—1)"-1! (m — 1)! (log 4)~?, a particular solution of (4) corres- 
ponding to A, = — m/2 log 3 — # is given by 
= ( ])™+r 
= —z\mtr 
yo = exp (—m/2 log + — 2)z Sipe Tin bral es Hugi 
r=0 
that is, 
z/2 Q-mz/2 ~ ERED a —.———(2-)' wALI 
oc Goat a ree Prt) U(m+r +1) 
r=0 
Replacing r by —m + r, we obtain 
Yo = ene? 2m (2) 





Tir + 1) (—m+r +1) Gel 


(equation continued on p. 454) 
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co 
mz /2 dD 272) 
i eee »; re+)P(—axntr¢)h @ 
r=( ° 


Thus we see that for any value of n 


z nz . See —z)r 
D, (z) =z e#/* 2 “> Tether ) pa br 4) 
and 
co x. : 
D_, (z) = e-#/2 Qnz/2 > TOFD ere ey ser g KY 


are solutions of (4). 


3. THEOREMS ON BESSEL TypE FUNCTIONS OF THE First Kinp 


Ifn is an integer, it is observed from (7) and (11) that D_, (z) = (—1)" D, (z) 
and so the solutions D, (z) and D_,(z) are not linearly indopendent.’ The following 
theorem strengthens the assertion. 


Theorem 1\—The Wronskian of D, (z) and D_,(z) is given by 
W (D,, D_») = log 3 e~ ae 
Proor: If P(z) and Q (z) are two solutions of (4), it is not difficult to show 
that 
W(P,Q)=A es 
where 4 is a constant. 


The proof now follows from the fact that 


-z —nz ae Whale —- 
D, (z) = e-#/2 2-nz/2 Patpn +22) 


and 


] 
vie = —2/2 nzf2 ae ee ee ps 


Remark 1 : Since W (D,, D_,) vanishes only when 7 is an integer, it follows that 


D, (z) and D_,(z) are linearly independent solutions of (4) except when 7 is an integer. 


In the case when n = 


. -k m, where m is a Positive integer, a solution of (4) is 
given by 
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— p-z/2 9—mz/2 ~ Cl)" ; -2z\r 
Duy fz) e781 2-eFF Ep ly eter + 1) (2*)", 


r=0 


To obtain a second independent solution of (4) we should apply Case III (p. 434) 
of Das and Lahiri*. A second independent solution of (4) can then be taken as (— 1)” 


oy 
ae evaluated at A = A,, where 


06 
= (A — A,)cy + = Co Fr; (2=7)" 


and c, r, are as in (8), (9) and (10). Hence asecond independent solution of (4) is given 
by 


D,, (z) = zDm (z) — e-*'® 2™/? (log 4)-} eee (2-7) 


r=—0 


eon ce (—1l¥ {—A (m — 1) + H(r)+(m + 7)} 5 
Tee Pes) > Se pal ra Dalia roe 2 ote 


where H(p) = 1+4+... + I/p, H (0) = 0. 


Ifn = 0, then A, = A. = 4. One solution of (4) is 


Dy (2) = el? >a je ey 


and the other independent solution, in view of Case II (p. 432) of Das and Lahiri’, can 
be taken as 


Dy (2) = (1 + 2) Dy (2) — 2e#!* (log 4) a H(n) (2. 


rol 
Theorem 2 (Recurrrence relations)—(i) "D, (z) = 2-4/7 [Dy-; (2) + Days (2) 
(ii) 2D‘ (z) + D, (z) = log 1 2-2/2 [D, -1(z) — Das (2)]. 


Proor : Multiplying D, (z) by 2~"/? and differentiating with respect to z, we 
obtain 


a (2-"/2 D, (z)) = — 4. 2°)? D, (z) + log + e77?? 
Z 


(equation continued on p. 456) 
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=) t. n 


Seats MCE IEE EKG 





= — $2-™? D, (z) + log? 27 D,_, @). 
This gives 
D, (2)+ Flog 4 D, (2) = — 3 D, (z) + log $. 2-*” Dy; (2). 
that is, 
D‘ (z) +( 5 log 4 + 4) D, (z) — log 4.2-#/? D,_, (z) = 0. ..(14) 
Again, multiplying D,, (z) by 2™/? and differentiating with respect to z, we obtain 
& (2! D(z) = — } 2"1? Dy (2) + log 3 e-*P 


< (- lyr —z)\r 
2TOrGtreH 


= — 12"! D, (2) + log t e#” 


co 
(— 1) (2: 2)s+1 Lnz E (n=1)z/2 
STE FDOT Gt itstiy — —F"" Pa @) — log t 2mm 
mr Da 2h 
This gives 
D', (z) — n/2 log 3 D, (z) = — $D, (2) —log 4.2-#!* D,,, (2), 
that is, 


Di, (z) — (n/2 log + — 4) D, (z) + log } 2-4/? D,,, (z) = 0. -..(15) 
From (14) and (15) we obtain 


nD, (z) = 2-*2[D _, (z) + Dass (z)) 
and 


2D), (2) + D, (z) = log 4, 2-4/2 [D,-1 (z) — D+ (2). 
This proves the theorem. 


Theorem 3 (Generating function)— 
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exp {—2/2 + 2-#"* (t— 1/t)} = E D,(2).1%,0<8 < | t| <o. 
n=—0o 


PROOF: We have 





and 
— 2-7/7 /¢) = Eel Pag (2-7/2) | 
exp (— 2-#//f) = 1 — ie tes ae ie 
Then 
oe) 
exp {2-#"( f= +) a > £"(2-2/2)s ee (Qany 
4 n+ 1 
n=— 00 
ae (272 /2)4 
2! (n + 1) (n ra 2) a a 
So, 


exp{— 5 + 278 ( Tha SS fern ane 


> (rey \ 
T¢+)lradtr+h 


= Sp. te) 1. 


n=—CO 


3 


This proves the theorem. 


Theorem 4 (Integral representation)—For integral values of n 
D,, (z) = = | e~2/2 cos (né an Q(l-z) /2 sin 6) dé. 
0 


Proor : In view of Theorem 3 


l ~z ] 
D, (z) = ioe e-7/? e? 3 (t — re f-*"!'dt 
Cc 


where C is a closed contour enclosing the point = 0. Taking C to be the unit circle 
| ¢ | = 1, we obtain 
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nw 
l a ™ 2-3/2 24 sin @ 
D, (2) = “Fa | e7zl2 e-ind 6 z/ sin do 
-% 


e~7/ cos (n@ — 2%-2)!? sin 6) dé 


I 
a |- 
ou 


the imaginary part vanishes since the integrand is an odd function. This proves the 
theorem. 


4. Besse. Type FUNCTIONS OF THE SECOND KIND 
It is clear that the function M,(z) defined by 


Mint D,(z) cos nm — D_,(z) ...(16) 
$10 Nar 

is a solution of (4) for non-integral values of nm. For an integer m, we define 

M,(z) = lim M,(z). 

n—>m 
Applying L’Hospital’s theorem, we obtain 
_ 1 § dD,(z) » aD_{2) } 
Mn(z) = + {eo «i Ce eee (17) 








The above relation further gives 
M_,(Z) = (—1)™ M,(z). 


For integral values of n, it may be noted that D,(z) and D_,(z) converge uni- 
formly in n and this gives that M,(z) coverges uniformly in m and consequently the 
function M,(z) is a solution (4) for any value of n. It is observed that this solution 
resembles the classical Bessel’s function of the second kind. 


Theorem 5 (cf. Theorem 2)—(i) n M,(z) = 2-7" [M -1(2) — M,4,(z)]; 
(ii) 2M, (z) + M,(z) = log }. 2-42 [M,1(z) + Mn,,(z)). 
PRooF : Following the proof of Theorem 2, we obtain 


d 
7 (2772/2 D_,(z)) = — } Q-nz/2 D_,(z) — 2-(w+12)2 log + D_n41(z); 


d 
dz 2%" D_y{2)) = 4 27? D_ (2) + 2-4/9 log 4 D_,_a(2), 


Utilising (16) and simplifying, we obtain, for non-integer n, 





d 2 
dz (2 "7" M,(z)) = — } 2-27 M,(z) + 2-™+1)/2 Jog 4 M,, -1(z); 
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< (2"2!2 M,(z)) = — $ 2”!* M,(z) — 20-072 log } Mys,(z). 


The rest of the proof is analogous to that of Theorem 2. We omit the details. This 
proves the theorem. 


We present now the series expansion of the function M,,(z) for integral values 
of m (> 0). 





m—) 
= ae ! 
Theorem 6—M,,, (z) oe + e~z/2 Qmz/2 SS (m - 1)! (2-2)° 
r=0 


“a _ enzl2 Q-mz/2 Geer Wiork — 8 1) (4 ED 


r=0 


where # (z) = I” (z)/T (z) and the first sum is equal to zero if m = 0. 


Proor : Since the series (12) for D,(z) converge uniformly in n, we can differ- 
entiate term by term with respect to ” and obtain 


co 





113 Oe SAA Ds Cay (y 
oe ica Pe ial Bae T(r + 1) l(m + r+]) 
r=0 
(= 1) yd 
ga il les ate 3 T(r + 1) al lim + r+ 1) } 
r=0 


Sol) a a 


Z —z/2 —mz/2 — 
.; z l0n.4,¢ are Tir + 1) Mim +r + 1) 
0 


r= 


oO 
—z/2 J-mz/2 oy 1)’ ty = 
ar T(r +1) Tim +r 42) 
r=0 





oo 


| 3 m+r+] \ 
x \1+ aaaT s(m+r+1+s) 


s™1 


lee) 
-2/2 Q-mz/2 Ss = ee 1)" (23 7 
ren Trt brm+rt)) 


r=0 





x {Flog} — vm +r + vt, 
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where y is the Euler’s constant. 


Similarly using (13), we obtain 





ee) 
dD_,(Z)_ / = e 7/2 Qnz/2 Ree en} wl 1 
dn ian Peo+1)%(-n4+r+41) 
ro 





x {— $ tog 4 tener ty tl 


Forr = 0, 1, ....m — 1, both (—n +r + 1) andb’(—n+,r + 1) tend to 
ec as n + m so that the first m terms of the last series become indeterminate. 


We note that 





a ws ne aryin oiee 
Mair De omnes 
so that 
w(-n+r+1) _ (-n+r+1) =. (n—r) sine (n—1r) 


P(=n-brpl) > (snr +e 
yx in —r) + xcot a(n — r) 
Hence for r= 0,1; ..; m— 1 


lim v(—n +r +1) 
n>m T(--n +r +1) = (— 1)"-"(m — r — 1)! . 


Thus 


aD_,{ 
cian ro as (S16 ee > apa ee ——___—_—— 27") 


co 


+ e7#/? Qmz/2 eee 
Pir+1)l(—m+rF 


fan 


x ee: 


= (— 1)" e713 zea S (m= Wt (i= M (g-2y. 


r=0 


—])™ e-2/2 mz (-—1) Z zjecte 
+ (—1)" e-#/2 2 6 Grae se 


(equation continued on Pp. 461) 
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x 1- 7 bed tart yb. 
Hence from (17), we obtain, for integral values of m (> 0) 


M w= 


M,(z) = — Ae en? /2 Qmz/2 S: a Soa ose (2-7) 
7 r! 


r=-0 


co 
l iaires (—1) (274) 
Bo gaze mz [2 
+ 7° 2 r!(m + r)! 


r=0 


x {elope — drt) —dimtr +t, 


where the first sum should be set equal to zero if m = 0. 


This proves the theorem. 


Utilising relation (12), we obtain an equivalent statement of Theorem 6 as 
follows : 


M,,(z) = ( x log 4 ) Z D(z) ae +. e7 2/2 Qmz/2 Slmiesr t (2-97 


r™=0 


(oe) 
Be eee oo vg Caly GY 
OS z/2 mz /2 BED one Me, ved eS de 
coh aie ritn +r)! 


r=0 


{b(r + + b(m +r + 1). 


5. REcIPROCITY RELATION 
Let D,(z) be defined by the formula 


exp (— 5). exp ( 2? ( 1—+))- -> D,(z)t",0 <8 @ | t| <ee, 


ratis) 


and let us call the function on the left, the generating function for D,(z). By actual 
computations we see that 


3 D,(z)t” = SS {exp (- 5) 7 ba el be Ss wae 
a™=—OO n=— OO r=0 


valid for all values of z and 0<5<|ft|< ee. This gives the Dirichlet series 
expansion of D,(z) as 
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Z nz/2 i 1)’ (2-7) ——— 
D,(z) — Xp ( — 3) genes, To+)ra¢r ae 1) ° ...(19) 


In view of Theorem 2, we have 


nD,(z) = 2-*!? [D,_,(2) + Dagi(z)]; -- (20) 
2 Di (2) + D,(2) = log $ 27? [D,-.(2) — Dys,(2)). .. (21) 
Eliminating D,(z) from (20) and (21), we obtain 
D‘ (z) + ( 3 log } + 4 ) De) = log } 2-*/? D,_,(z). .. (22) 
Differentiating (22) with respect to z 
Dy(2) + (4 log + 4) D,() = 4 (log §)* 2-2” Dy_x(2) 
+ log $. 2-#2 D’_(z), 


Substituting for D)-, (2) from (22), (n replaced by n — 1), we obtain 


D, (2) + (5 logs +4) D,G)= 4 (log dy 21 Dy-a() 
+ Glog 4) 2° Dy-s (2) 
~ tog 42-1" ( "5" toga + 4 ) 
x D,- (2). 


Eliminating D,-»(z) from the above relation and (22) with n replaced by n — 1, we get 


Date) + (5 logs +4) Di) = (wm — 1) (og YY2” D,. 





— (log 4)? 2-* D,(z) — log § 2-712 ( aa ? tog 3 + +) D,-; (2) 


= (n — 1) log 3 | D‘ (z) + ( yloet +4 )D, «} 





n— 2 
~ (734 24+2) ne - ( Flos t +4) 


a ss 
x ( 5 log} +2 — log + ) D,(z) — (log 4)? 2- D,{z), using (22). 
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This gives 
Ds (z) + D(z) + | (log 2) r+ —(F logt + 4) (4 tog} - } 
x. D(z) = 0. ules) 
This shows that D,(z) is a solution of the equation 
ae ra ie +, | = Rees rar Be » wen (24) 


where « = log + and n is an arbitrary constant. 
The differential equation (24) is eqn. (4) with which we started. 


Remark 2 : If in the equation (4) we take « = log 1/p where p is a positive 
integer (> 2) instead of log 4 then except longer calculations, analogous results would 
be obtained. 


REFERENCES 


1. L. Cesari, Functional Analysis and periodic solutions of nonlinear differential equations. 
Contribution to Differential Equations, John Wiley, 1 (1963), 149-87. 

2. L. Cesari, Mich. Math. J. 11 (1964), 385-14. 

3. A. G. Das, and B. K. Lahiri, Rend. del Circolo Mat. di Palermo. II, Tomo XXXIIL (1984), 
425-35. 


Indian J. pure appl. Math,, 19 (5): 464-476, May 1988 


GROWTH AND APPROXIMATION OF GENERALIZED BI-AXIALLY 
SYMMETRIC POTENTIALS 


G. P. KAPOOR 
Department of Mathematics, Indian Institute of Technology, Kanpur 208016 


AND 


A. NAUTIYAL 


Systems Analysis Division, Defence Research and Development Laboratory 
Kanchanbagh, Hyderabad (A.P.) 


(Received 22 June 1987) 


The present paper deals with growth and approximation of solutions (not 
necessarily entire) of certain elliptic partial differential equations. These 
solutions are called Generalized Bi-axially Symmetric Potentials (GBSP’s). 
The GBSP’s are taken to be regular in a finite hyperball and influence of the 
growth of their maximum moduli on the rate of decay of their approximation 
errors in both sup norm and L®-norm, 1 < 8 < ©, is studied. 


1, INTRODUCTION 


Generalized bi-axially symmetric potentials (GBSP’s) are the solutions of the 
elliptic partial differential equation 


o°7H o°H 2a+1 dH te 20+1 0H 


“ae + Be P Oy Pe age Oi ett 
wt Let) 
which are even in x and y (Gilbert*). A polynomial of degree n which is even in x and 


y is said to be a GBSP polynomial of degree n if it satisfies (1.1). A GBSP H, regular 
about origin, can be expanded as 





co 


H =H (r,9)= Sy a, r* pst (cos 29) . ent Bea) 


n=0 


where, x = r cos 6, y = rsin é’and ee (t) are Jacobi polynomials. 


Let Dr = {(x,y): x? + y? < RY, 
A GBSP H is said to be regular in Dp 
compact subsets of Dr. 
R' < R but for no R’ > 
Recently, McCoy 
GBSP polynomials and 


0 < R<coand Dg be the closure of Dr. 
if the series (1.2) converges uniformly on 
Let Hr be the class of all GBSP’s H regular in Dp, for every 
> RK. The functions in the class H. are called entire GBSP’s. 
”8 considered the approximation of an entire GBSP H by 
found the rate of decay of approximation errors 

E, (H, 1) = inf || — gi, = inf { 


a max | H (x, y) |} 
ger, 7). 


Er, (se y) Ee D, 
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and 


E,,8 (H,1) = inf ||H — glh,s 
gEty 


= inf {JJ u(x, y) | (x,y) | & dx dypls 
gGrn D, 


» a weight function, 1 6 < oo, in terms of the growth parameters associated with 
the maximum modulus function 


M (r, H) = max | H(r, 6) |. 
2] 


The effect of rate of growth of M (r, H) on the coefficients a, in (1.2) of an entire 
GBSP H were studied by Fryant®. 


However, for a GBSP H which is not entire, no attempt seems to have been 
made so far to study the approximation errors or the coefficients in (1.2), vis-a-vis, the 
rate of growth of M(r, H). The present paper is an effort in this direction. 

A GBSP H is said to be regular on Dz, , 0 < Ry < oe, the closure of Da, , if it 
is regular in Dr, for some R’ > Ry. We denote by Hr, the class of all GBSP’s H regu- 
lar on Dr,. For H € Hr,» S¢t 


| Hk = max |4H(x,y)| 
a Le De 


ow 


( Jw (Ro, 8) | H (Ro, 8) | ® dé yn iss Bieese tI) 


0 


\| H las 


I 


([f emriweni tard J 1ctce 09 


Dr, 


|| H las 


where the functions w and # are positive and integrable (in the sense of Lebesgue) 


such that 1/w and 1/@ are bounded. Then || - |lx is the uniform norm on Ha, whlie 
0 


ll: Ie vg and Il - Ili _s are L*-norms on Hr,. For HE Hr, the approximation errors 
0 0 


E, (H, Ro), E.., (H, Ro) and E,,, (H,Ro) are defined as 


E,(H,R,) = inf | H—8 lz, tL) 
gE 
E\,, (H,R.) = inf ||H—8 eet (1.6) 


gEtn 
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E.,, (H, Ro) = inf |H— giles oi E78. 


gern 
where 7,, consists of all GBSP polynomials of degree at most 2n. 
A GBSP H € Hr, 0 < R < 9, is said to be of order P if 
log* logt M (r, H) 
= p (H, R) =! et 
Pee PE RE ee 


where, log+ x = max (0, log x). If0 < Pp < co, then the type 7 of H is defined as 


ae a logt M (r, H) 
T =T(H, R) = Heusen “(RAR —r)yP 


In Section 2, we obtain characterizations of rate of decay of approximation 
error E, (H, Ry) interms of order P and type T of H € Hr, 0 < Ry < R <oo, The 
Characterizations of rate of decrease of Re. (H, Ro), i = 1, 2, and a, in terms Of these 


growth parameters of H € Hp are found in Section 3. 


AEE 


2. APPROXIMATION Error E, (H, Ro) 
We need the following lemmas : 


Lemma |1—Let H € Hr, R > Ry. Then, there exist GBSP polynomials g, € 7, 
such that 


|| 1 — 8, liag < KM (r, H) (n + 1)e*?? (Ry ir)2e+)) 
for all r sufficiently near to R and all sufficiently large values of n. Here K is a constant 
independent of m and r and g = max (<,). 


Proor : Set 


£, = SS a, r2* P™®) (cog 24) , 
k=0 


Then g, € 7,. Using (1.2), we get 


co 
|| H—g, llr, < >» | a, | R°* | Pk”) (cog 24) | 


k-n+1 


| oo 
<Page D, eR ket e+ DPE +D 


k=n+) 


selon) 
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since, we have® (p. 168) 


(« 5B) C(k + 1) 
pe Cpe ctl Mr SPRINT TAT » q = max («,8). 


A222) 
For H € Hr, we have ({[3, (8)]) 


[(2k +a + B+ 1) A(k, a, B)A(a, BP ...(2.3) 


for every r < R, where 


—_ VT (kK+1) 0 (k+a+8+1) — F(a+i)r(p+l) 
As Sy Pd P(k+at+1) 0 (kK+8-+1) ’ Ries P) coy I (a+6+ 2) 
2.4) 


Combining (2.1) and (2.3) we get 


M (r, H) T(k+q+1) 
| H — galley S T(q+l) V A(z, Ta S mel (eaeiye 


ken+1 
x [2k +a+8 +1) A(k, a, B® (Rolr)2*. — ...(2.5) 
Since [ (x + a)/I' (x) ~ x* as x — 9, we have 


lr (k+q+ 1) 
T (k+1) 


as k — cc and so 


I (k+q+1) 
T (k+1) 


for allk > ky. Thus, forn > ky andr > r*, wherer* = (R + Ro)/2 if R < ee and 
= 2R, if R = oe, using (2.5) and the above inequality, we get 


((2k+a+B+1) A (k, a, 8)? ~ V2 ke? 


((2k+a+B+1) A (k, a, BN? < W2W2keHl 


co 


M (r, H) YP Par Ve ka+1/2 (Rolr)2* 
| 7 — gnilz, < Tepe ie 2 V 2A (a, 8) >, (Rolr; 
M (r, A) WFR Ve +1/2 2(n+1) 
— = 2 2A (2, + 1)2*4/? (Rolr) 
< (q+) 2 2A (a, B) (n 0 


oo 


41/2 
x SS ( 1 + Ear) (R,/r*)*. 


k=0 


The lemma now follows from the above inequality. 
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Lemma 2—Let H € Hr, . Then, forn > 1, we have 


| a, | Ro" < 2 ((2n+a+8+1) A (a, B) A (n, @, B))/? E,-, (H, Ro) 


where A (n, «, 8) and A (a, 8) are given by (2.4). 


Proor : From the orthogonality of Jacobi polynomials® (p. 68) and the uniform 
convergence of the series (1.2) on Dr, we have 


a, Ri" | ((2n+a+8+1) A (n, «, B)) 


T/2 
= 2 [ H (Ro, 0) pe (cos 26) sin?**? 9 cos*®*! g dg. 
0 


Thus, for any g € 7,-, we have 


a, R®" / (2Qn4+a+8+1) A (n, , 8)) 


W/i/2 


= 2 | (H (Ro, 8) — g (Ro, 8)) peat (cos 29) sin***! 9 cos28*1 9 do. 
0 


...(2.6) 


Using Schwartz’s inequality and the orthogonality of Jacobi polynomials, the above 
relation gives 


| an | Ry” < || ~glle, (2n+0+8+1) A (n, «, 8) A(z, Ba, .(2.7) 


By the definition of E, (H, Ro) there exists a GBSP polynomial g € 7z,-, such that 


| 1 — g lla, < 2E,_, (H, Ro). ; --. (2.8) 


Taking, in particular g = & in (2.7) and then using (2.8), we obtain the lemma from 
eae 


Using Lemmas | and 2 we now prove : 


Theorem 1—Let H € Hp, . Then H € Hr, R> Ro, if and only if, 


lim sup (E, (H, Ry)!" = RR. 
ll co 
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PROOF : First, suppose that H € Hr. Then using (1.5) and Lemma 1 we get 
E, (H, Ro) < ||H—g,|k, <K M (r, H) (n + 1)9*1/2 (Rofr)2+), ...(2.9) 
The above relation gives that 


lim sup (E, (H, Ro))!/?" < R,/r 


1c 


for all r sufficiently near to R and so 


lim sup (E, (H, Ro))'2" < Ro|R. 
71-—> 00 


On the other hand, using (2.2) and Lemma 2 we get 


oo 
> a, rn” pe (cos 26) | “= 
n=O 


FA eee Cn YAS PS 2.8 RD ( Fe) 


Pr (n+q+1) 


ci P (n+) P@+) 














Pr (q+!) 
x (Qnta+e+1) A (n, a, Eye SE EIED (2.10) 


It follows from (2.10) that if lim sup (£, (H, R,))'*" < R,/R then the series (1.2) 
n> co 
converges uniformly on compact subsets of Dr, for some R’ > R and so, H is regular 


in Dd,’ which is a contradiction. Hence, lim sup (E£, (H, R,))'/*" = R,/R. This 


n-->oo 
proves the necessity part of the problem. 


Sufficiency part can also be proved similarly. This proves the theorem. 


For proving our next theorems we need the concepts of order and type of a 
function of a single complex variable u analytic in the disc |u| < R. 


Let f(u) be analytic in | u| < R, 0 < R< cc. Then the order Py of f(u) is 
defined as 


log* logt m (r, f) 
Po = lim suP ~Tog (R(R—P) 


where 


m(r, f) = max | f(u) | 
}u|<r 
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and if 0 < P, < oo then the type 7, of f (u) is defined as 


— togt mrs f) 
(ND SUD CR ree 


We now have : 


acd . . 
Lemma 3 (MacLane’, p. 47; Beuermann’)—Let f(u) = z b, u" be analytic in 
Ju |.< R(O < R < co) and has order Py (0 < Py < ©). Then, 


Po = lim sup pee ah a Ee 
: aon. 680k — log” log p ike 


9° . . 
Lemma 4 (Kapoor’, p. 256; Bajpai et al/.1)—Let f(u) = ve b, u" be analytic in 
|u| << R(0< R< oo), Then/fis of order Po (0 < Py < 00) and type 7», if and 
only if, 


‘73 
Vo = Pork 1) oo Is 
Po 
Py 
where 


+ b R ork 
v) = lim sup Mogl Via) a.. = ) 


n-> oo n 0 


satisfies 0 < vy < oo, 
We now prove the main theorems of this section. 


Theorem 2—Let HE Hr, 0<Re< co, be of order P and Ry < R. Then 


: log* log* (E, (H, Ro) (R/-Ro)*") 
a} ; ; Sok eee | 
f= Thm SUP “Tog n — log* log* (E, (H, Ro) (RIR,)*) (3) 


Proor : Let the limit superior on the right hand side of (2.11) be denoted by d. 
Obviously 0 <d€ co, First let 0 < d < co and d’ be an arbitrary number such 
<a 


thatO < d' <d. Then, by the definition of d there exists a sequence {n,} of positive 
integers tending to co such that 


an. q itd’) 
log (En (H, Ro) (R/Ro) ") > nf (2.12) 


forkrea ae Using (2.9) and (2.12) we get 
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log M(r, H) > ni !@*") 4. Ing log (r/R) — (q+1)]2) log (m+1) — log K 
2.(2:13) 


for all sufficiently large values of k and all r sufficiently near to R. Let {rx} be a sequence 
defined by 


ny = {3 log (R/r,)}- O74” aa(2.1A) 


then r, ~ Ras k + oe, Now, using (2.13) and (2.14), for all sufficiently large values 
of k, we have 


log M (ra, H) > ~areqr (log (Rira))-#” + (q-+1)2) (1 +4’) log (3 log (Rira)) 


— log K+ 0()). 
Since log (R/(R — rx)) ~ — log log (Rir,) as k + ce, the above relation gives that 


é log* log* M (rx, H) ~ 
lim sup —~>—>-———.. © d’. 
See Tog (RI(R — rx) 
Since d’ (< d) is arbitrary, this in turn gives that 
P > d. vos(eul DD 


Obviously, (2.15) holds for d= 0. For d = co the above arguments give that 


Pp = co, 


On the other hand, using (2.10) we get 


2A («, B) 


M (r, H) < | a) | + Path) m (r, h) ..(2.16) 
where, by Theorem 1, 
< UD 1 
hw) = > (Qnt2t B41) A Cn, «, sys ETE? 
ny 
x E,-1 (H, Ro) (ul Ro)" eek a) 


is analytic in |u| < R. Using (2.16) and applying Lemma 3 to h(u) we get 
P <d. Hee YY 
Theorem 2—now follows from (2.15) and (2.18). 


Theorem 3—Let H € Hr, 0< R< 2&9, and Ryo < R. Then A is of order 
p (0 < P < oo) and type T, if and only if 


y = (p+ 1) (2/0) T 
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where 


(log* (E,, (H, Ro) (R/R,)*"))e*? 


v = lim sup per 


n> co 


pe SEN 


satisfies 0 < y < co, 


Proor: (i) First, let H € He be of order P and type T T(< ©), Then, for 
e > 0, (1.8) gives that there exists rp = ro (€) such that 


log M(r, H) < (T + €) (RI(R — r))e ...(2.20) 
forro <r< R. Using (2.9) and (2.20) we get 
log* (Ey (H, Ro) (R/Ro)"") < (T + €) (R/R — r))e + 2n log (Rr) 
+ (q+1/2) log (n+ 1) + logt K ae ty iy 


for all r sufficiently near to R and all sufficiently large values of n. Choose a sequence 
{r/ } as 


RR ~ 7) = (5 eee (2.22) 


Clearly r, + Rasn-— cc, Using (2.21) and (2.22) we get 


T+e)(1+ p) (2n)e/(1+e) 
log* (E, (H, Ro) (RIR,)*") < ee (1+P+0(1)) 
for all sufficiently large values of n. 


On proceeding to limits, the above inequality 
gives that 


v < (P+1)F+1 (2/p)e 7. 223) 


y using (2.16) and applying Lemma 4 to the 
This proves the necessity part of the theorem. 


The reverse inequality in (2.23) follows b 
function h(u) given by (2.17). 


(ii) Suppose that 0 < y < oo, Then (2.19) gives that 


; log* logt (E,, (H, Ry) (R/R,)*") 
p = ] et Pr E BASSES) Fi, **0 “a ee: 
noc) log n—log* logt (E, (H, Ry) (RIR ™) 


and so, by Theorem 2, H is of order P. Sufficienc 


y part now follows from the necessity 
part of the theorem. 


This proves the theorem. 


l 9 
3. THE Errors Eg (Eo a). Ews (H, Ro) AND THE COEFFICIENTS a, 
In this section we obtain some results analo 


gous to Theorems 1, 2 and 3 for the 
approximation errors in L&- 


norms and the coefficients a, in (1.2). Since the techniques 


a 
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used in proving these results are the same as those used in Section 2 we only give the 
outlines of the proofs. 


Let H € Hr and eixeR: Using Lemma 1, (1.3), (1.4), (1.6) and (1.7) we get 
Ee (A, Ry) << Ki (n + ])e+1/2 (R,/r)2rD M (r, rhe i= he PS 2 KT 


for all r sufficiently near to R and all sufficiently large values of n. Here K, is a 
constant depending on R), w and 8 only and K, a constant depending on Ry, # and 5, 


On the other hand, by (1.6), for HE Hr, there exists a GBSP polynomial 


g* , © 7,_1 such that 


1 (eee (7, Ry) = || i g* 


1 
Shp & 
0 


ow 
> pie ( [170 — et, (oa) | Fao) 3.2 


since I/w is bounded we have w aA/t, T>0, For 8 >.1 shoose 7 such that 
(1/5) + (1/n) = 1. Using Holder’s inequality we get 


on 
| | H (Ry, 6) — * , (Ry, 8) | dg 
0 


QW 


<( Jie 0 - et, coo ia) ( [aon 233) 


0 


Combining (3.2) and (3.3) we get 


27 


| | H (Ro, 0) — 25 (Ro, 6) | de 


0 


1 


2En-vs (H, Ro) > T#I8 (27) 


| 1 HR 0) — ef (Ry dd. 3.4) 


0 


4 


for 5 > 1, since GBSP’s H and g,_,areeveninx andy. For 3 = 1, (3.4) is obvious 


with y = co, Taking, in particular, g = g* , in (2.6) and using (2 2) we get 
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| a, | Ri" | ((2n + « + B + 1) A(X, ag, B)) 





1/2 
ae eh 
0 lace) 
Combining (3.4) and (3,5) we obtain 
lan | RK T? ® (27)!/" (2n+a4+8+ 1) A (n, a, B) PP (nt+gt+!) 
c a lr (q+1) [ (n+1) 
xX pha (H, Ro). ae (3.6) 


Using (3.1) with i = 1 and (3.6) we see that E,(H, Ro) can be replaced by 
E,,, (H, Ro) in Theorems 1, 2 and 3. 


Similarly, by (1.7), for H © Hg,, there exists Z,_; © 7,-1 such that 


2E je (H, Ry) > WH — Bn—alle, ; 








J l 1/8 
ae (|| H (X,Y) — Bn-1 (x, y)|® dx dy) 
Dr 
0 
—_ l 
x JJ | (x,y) — Baa (x, y) | dx dy en) 


Dr 
0 


where > iy, T= 0 and (1/8) + (/yy}= 1. . From the orthogonality of Jacobi 
polynomials and uniform convergence of the series (1.2) on Dr, we have 


ayr'"((2n+a+8+1) A (n, «, B)) 
t/9 


=— 2 | (H (r, 6) — Bn-1 (r, @)) PS” (cos 26) sin2*? 9 cos?8t! 6 do 


0 
forO<r< Ro. Using (2.2) we get 


ly, [:r@h (nee B+1) A (n, «, B)) 


— TP (n+q+1) 
< 2 @+l) P+ | | H(r, 0) — &,-1 (r, 0) | do 
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since H and g,-, are even in x and y. Multiplying both sides of the above inequality 
by r dr and integrating from 0 to Ry we get 


| a, | Ro"*? [((2n+2) (2n+a+841) A (n, «, 8)) 


T (n+q+1) _ 
<arG@t+h Path {| | H (x,y) — Bn-1 (x,y) | dx dy. ...(3.8) 


Dp 
0 
Combining (3.7) and (3.8) we obtain 


la, | Ret? Poe Bae (7 RG)" (2n+2) (2n+a+B6+1) A (n,a,8) P (ntgt+) 
Stemi Pr (q+) T (n+1) 


x Sale (H, R»). my 


Using (3.1) with i = 2 and (3.9) we see that E, (H, R,) can be replaced by 
E;,, (H, Ro) in Theorems 1, 2 and 3. 


Finally, let H € Hr. Then, forr < R, using (2.10) we get 


M(r, H)<|a.| +m(r, AIT (q+1) ...(3.10) 


where, by Theorem | of Fryant® 


T (n+qg+1) 


h(u) = T(nt+l) 


| a, | wen 


n=1 
is analytic in |u| < R. 
Using (2.3) and (3.10) we obtain the following coefficient characterizations of 
order and type of H © Hr, 9< R< oo: 
Theorem 4—Let H € Hr, 0 < R < ©°, be of order p. Then 


= __log* log* | an | R” _ 
p= tim SUP jog # — log* log* | a, | R™ 


Theorem 5—Let H € Hr, 0< R< oe. Then, His of order P (0 < P < ©0) 
and type 7, if and only if 


Ve = (P+ I)e*? (2/p)e T 


where 


(log | ay |_R*)F* 


v, = lim sup a 


n>o 


satisfies O < vy < ©. 
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In this paper, the author has investigated sufficient condition, which cannot be 
relaxed to any appreciable extent, for absolute summability factors for | R, 
exp {n/log (n + 1)},1|,for both general infinite series and Fourier series, 
which consequently improves the corresponding results due to Lal [Proc. Am. 
Math. Soc. 14 (1963), [311-19]. The class of functions, considered by the pre- 
sent author, for the Fourier series strictly contains the class of functions con- 
sidered by Lal.® 


1. DEFINITIONS AND NOTATIONS 


Let f be 27-periodic and L-integrable over (— 7, 7) and let, without loss of 
generality, the Fourier series of f at x be 


co io 2) 
Dy (a, cos nx + b, sinnx) = Ss Ae 2): 
n=l n=) 


Let x be any fixed real number. Then we write 


dt) = 41f (x +O + f(x — 0} (1.33 
p(t) = t-? f ¢ (u) du ry ee 
9 
P(t) = $(t) — oi (t) atlia 
Pr (1) = $s (0) — 8 f dy) du (4) 
e(n) == exp {n/log (n + ID} (n > 1) eat Ls) 
p= f-* Yalog (2 + 1) --+(1.6) 
Az, = Zn — 2n41 (for any sequence (Zz), Cd @ fe a 


We also write 7 = [27/t], the greatest integer less than or equal to (2z/t). 


It is known that the summability methods | R, e(7), 1 | and | R’, em), | | are 
equivalent (See Mohanty’, footnote to the page 298). 
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A series & d, € | R, e(n), | | if and only if (see Chandra’) 
n=1 


> Adle@ | Se (m dy | <x. (18) 


n=1 m=1 


Let L (—n, 7) be the set of all 2n-periodic and L-integrable function f. Then we 
write 


S(f; di) = {f © L(—x, 7): 4, (t) © BV (0, m)}. 
Given a function 4, we also write 

GU Brea 
for the set of all functions f € L (—n, 7) for which 

lim, 

y>0+ J ¢: (u) du (1,9) 
exists (in the sense of Cauchy) and 

(P: (t)/h (t)) €© BV (0, 7). .-.(1.10) 

In the case h(t~!) = 1 for all t, we simply write C ( f; P;) for C ( f; P,; 1). 


2. INTRODUCTION 
Lal® proved the following : 


Theorem A—Let the sequence (y,) be 


convex, and eat och) 

Eon ea kee 33-3} 
Then 

SES hws ss Ay (x)¢, € | N, Un + 1)]. ...(2.3) 


In an attempt to prove Theorem A, Lal® also proved the following : 


Theorem B—Let (y,) satisfy (2.1) and (2.2) and let 


> in 
Son en (2.4) 
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Then 
co 
ud en E LN, TQ -F:1). |. 


We first observe that if (y,) satisfies (2.1) and (2.2) then 


O<Ay, J}, -+.(2.5) 
2) 
% log (n+ 1)Ay, < ©. of Ci) 


For (2.5), see Zygmund"’; p. 93 and Chow’ and for (2.6), see Mohapatra!®. Also 
from Lemma | of the present paper it follows that the summability method | R, e(n), 1 | 
is weaker than | N, 1/(n + 1) |. Thus, in view of (2.5) and (2.6) which are conse- 
quences of (2.1) and (2.2), we replace (2.1) by 


co 
x n-1 | Ay, | log | (n + 1) < © an f2cs) 
and the method | N, !/(n + 1) | by | R, e(n), | | to improve Theorem B in the follow- 


ing form : 


Theorem 1—Let (y,) be a given sequence satisfying (2.7). In order that (y,) 
should be such that, for every (d,) satisfying (2.4), 


Sees uence (sists (2.8) 
Aa 


it is sufficient that (2.2) should hold. 


Remark 1: The following example shows that (2.2) cannot be relaxed to any 
appreciable extent : 


Example—Let 
Yn, = I/log (n + 1) and d, = (-1)"(" > 1). 
Then (2.2) does not hold but (2.7) holds. Also 
> mdm = >> (—1)"m = O(n), 
m=) m=1 


which proves that (2.4) holds. Now to prove that (2.8) does not hold, it is sufficient to 
show that, in view of (3.1), 


ioe) 
x de, &©|N, Win +1) I. 
n=1 
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Suppose, on the contrary, that 


S de md (1 n € [Nie 
n=1 


n=) 


then we must have (see Corollary to Theorem | of Das®) 


> 


n=) 


(hin e 
n log (n + 1) 





< o, 





which is not true. Therefore (2.8) does not hold. 


This proves that if you relax the condition (2.2) to some appreciable extent then 
there exist sequences (d,) and (y,) satisfying (2.4) and (2.7) respectively for which (2.8) 
does not hold. 


Now, in view of Lemmas | and 2, we improve Theorem A, by proving 


Theorem 2—Let (y,) be a given sequence satisfying (2.7). Then 
[oe] 
SEc (LP?) >= 2 A, (x)e, € | R,e(n), 1], .. (2.9) 


provided (2.2) holds. 


This may be deduced from the following general theorem. However, a simpler 


proof of Theorem 2 may be obtained by combining the arguments of Remark 2 (given 
after Theorem 3) and Theorem 1. 


Theorem 3—Let (),) be a given sequence satisfying (2.7) and let there exist a 
function h (1/t) > 0 such that, uniformly in 0 < t < 7, 


“wh(i]o) ~ Oth(II)} va.(2.10) 
t h(1]t) log (n/t) = O (1) , Eve hth 
A(I)E wm] yi =O(1) (2.12) 
h (1/t) = n*| Ay, | log(n + 1) = O(1) E18) 
h (ye) n | Ay, | e(n) log(n + 1) = O {e(T)}. v.(2.14) 


We further suppose that as ¢ increases 
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d d 
07a, 2% | t ee (th cry} is non-decreasing for some 8 with 


0O<B<! w..(2.15) 
hfe HH t a @hauny} is non-increasing ..-(2.16) 

and 
Kal ; + (tA (1/t)) \] = 0 {h(n)}. (2.17) 


t=l/n 


Then, since (2.2) is included in (2.12) we have 
ie 2) 
fEC(S; Py h) => 2 A, (x) «, € | R,e(n), 1|. ...(2.18) 


Remark 2: Suppose h is such that C (f; P;) C C (f; P;; 4) then condition (2.2) 
cannot be relaxed to any appreciable extent. To prove it, we have as in Theorem 2 of 
Chandra and Gupta’. 


T 


A, (x)=2P, (m) cos nz +2 | P, (t) 1S Soe — cos nt) dt. 
0 


Now, we take f € C(f; P,), then f € C (/; P,; A). However, for f € COP); 
7 


A, (x) = O(n") + —| t cos nt dP (t) 


0 


and hence 


Sm 4.) = O(n) + 2 t dP, © (> m cos mt) 


m=) 0 m=1 
= O(n), 


which proves that (2.4) holds with dn = A,,(x). Thus by taking d, = A, (x) in the 
example of Remark |, if follows that (2.18) does not hold. 


3. LEMMAS 
We require the following lemmas for the proof of the theorems : 


Lemma |1—The inclusion 


| N, 1/2 + 1)| D | R, exp {nflog (n + 1}, 1 | TAC AD 
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holds. However, for 0 < B < 1, the inclusion 
| N, 1f(n + 1)| C | R, exp {n/(log(n + 5))®}, 1 | =s(3.2) 
is false. 


For the proof of (3.1), see Dikshit? and for (3.2), see Bosanquet and Das’: 
Theorem J. 


Lemma 2—The inclusion 
S(f; ¢1) C CCF; Pi) sh) 
is strict. 


ProoF: The inclusion (3.3) is simple to get. Thus we prove that the inclusion is 
strict. To prove it we consider an even function f € L (—71, nm) so that at x = 0, 


H(t) = f(t). 


Now we define 


d 
ne \ir (tlog (2n/1)) (0<t<7) G4) 


0 (t= 0). 
Then 
$1 (t) = log (2n/t) 
and hence f & S(f; ¢1). However, for f defined by (3.4), we obtain that 
P,(t) = —1., 
This proves that f € C (/; P:) and hence the lemma follows. 


Lemma 3—Let that function A(1/t) satisfy (2.15) through (2.17). Then, uni- 
formly inO0<t<vz, 


sin nu d d ; 
[Sone a u 7 (uh (ilu) du = O {n-* h (n)}. 


This may be obtained by using the technique of Lemma 3 of Chandra and 
Mohapatra‘. 


4. PRoor OF THEOREM 1 


co 
The series 2 d,e, € | Rie (n), 1 | if 
n=) 


y ioe) oo 
= 2 A (Ie (n)) | 2 (m) dm mM ¥mlog(m + 1) | < , 
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Now, by Abel’s transformation 


n n~l m 
= e(m)dmm y,, log (m+ 1) = = A Ym 2 k~ e (k) dy, log (k + 1) 
m=) engl | 


m1 


ys ae. k-! e (k) dx log (k+1)...(4.1) 


And the sequence (n~* e (m) log (n + 1)) is ultimately increasing with n, therefore by 
Abel’s lemma 


|£ mz e(m) d, log(m + 1)| = O {ne (n) log (n + 1} (4.2) 


m=) 


whenever (2.4) holds. Hence by using (4.1) and (4.2) in = we get 


= = O (1) = A (l/e (n)) 2 A y,, | m-} e (m) log (m + 1) 


+ O(1) = A (l/e (n)) n“ e (n) y, (») log (n + 1) 


= 4, + 2,, say. 


By changing the order of summation, we get 
2) 
Xi = O(1)=2 | Ay, | m log(m+ I)<-, 
pans. 


by (2.7), and using the inequality 
| Ae(n) | log (n + 1) = O fe(n + 1}. 
we obtain that 
co 
ret tl) Sn Pe | < ©; 
n=) 


by (2.2). 


This completes the proof of Theorem |. 


5. PROOF OF THEOREM 3 
Proceeding as in Theorem 2 of Chandra and Gupta*, we get 
T 


zi 
An(a) — &-[te, com apa cm) G-(SSee — co me) a 





+ 2 P, (x) cos nz 
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~ = 2 fie Faas (func J 


0 


- J fare, (t)/h (1/t)} [ A(1/t) (= — cos nt ) 


t 
sinnt d sin nu d d 
+ nt S(t A(t) — pa a | u <(u h(1|u)) } fe } 
0 
since, by (2.10), (2.11) and (2.15), we obtain that for 0 < 1 





(5-4) 
‘ ’ 
[u A(\/u) 5 (a — cos nu ) au 
0 


t 


= t h(1]t) (2 sine S — cos nt) oe \a(= ons 
0 


— ee, 


= jus u — ye A(1/u)) du 





= t h(1/t) ( sant — Cos nt) + sae 


d 
Senate bt) 


‘ 
__ [ sinnu d 
aE. - ian (u A(1/u)) I at 
0 
Now, by (2.18) id | {Pi (t)[h (t-)} | < 0 therefore by Lemma 3. we 


get 


A, (x) = O {n> h(n)} 4 ih (Pi (Ih (1]0)} th (Ife) cos nt 


0 


= [a (Ps (Ih (je) {ram+ Lon cep} si 


= A” (x) 4+ AY (x) + Ae (x), say 
However, by (2.11). 


S| A (x) €, 


n= 





= 00 Setiyy, 


n=) 
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co 


which is finite by (2.2). Thus 


at (x) €, is absolutely convergent. Hence by the 
absolute regularity of | R, e (n), 1 | 


co 


Ss A® (x) e,& | R,e(n), 1]. 


n=) 


Thus, whenever (2.10) holds and f € C (f; P,; h), for the proof of the sufficiency part, 
it would be sufficient to show that 


S A (l/e (n)) | y e(m) m™ ym log (m+1) K (m, ft) | 
n=l m=) 


= ({t A(Q1/t)}-*), (5.2) 


F . ee sin mt 
uniformly inO < t < 7z., where K (m, ft) is either cos mt or ee 


co 
We split up the sum 2 into & and & . Then by using | K (m, t) | < 1 and chang- 
a=) n<T n>T 


ing the order of summation, we get 


et Shayne) y, log (m+1) = O {(t AC/)}, .. (5.3) 


ne m=) 
by (2.2) and (2.11), uniformly inO <t<7. 


Now, by Abel’s transformation and Abel’s lemma, 


|S e(m) m™ ym log (m+1) K (m, t) | 
m1 
<E |Ay, Ji e(r)r™ log(r+1) K (7, t) | 
mre r=) 
+ | yn || > e(r) r-} log (r+1) K(r, t) | 
r=) 


= O(t") >> | Ay, | e(m) m- log (m+ 1) 
m= 
+ O(t 3) | yn | e (n)n7 log (1+ 1), 
since (n-} e(n) log (n + 1)) ultimately increases withn and 


$ Kin.) =O0UNO<a<b<~). 
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Hence 
ES =O(ift)E Ale (n)) = | Aym | e(m) m— log (m+1) 
a>T n=7 m=) 
+ O(I)t)E A (le (n)) | yn | e (2) n-} log (n+1) (5.4) 
However 


e(n) log (n+ 1) A (1/e (n)) = O (1), 


therefore the second term on the right of (5.4) does not exceed 
O()E nm | y, | = O (ACO), (5.5) 
n=T 
by (2.12) and the first term does not exceed 


O(1/t) > A (l/e (n)) y | A Ym | m-} e (m) log (m+1) 
nov m=1 


+ O(I/t) E A (le (n)) = | ¥%— | mt e (m) log (m-+1) 
n= m= 


O{(t e(T)* EA yg | mo? e(m) log (m+1)} 


+ O(1/t)E | Ay, | m7 log (m+1) 
m=T 


O (th (A/t))}, ...(5.6) 
uniformly in 0 < ft < 7, by (2.14) and (2.13). 
Thus combining (5.3) through (5.6), we get (5.2). 


This completes the proof of Theorem 3. 
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UNSTEADY MIXED CONVECTION LAMINAR BOUNDARY—LAYER 
FLOW OVER A VERTICAL PLATE IN MICROPOLAR FLUIDS 
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The combined effect of forced and free convecticn on the unsteady laminar 
incompressible boundary-layer flow of a thermo-micropolar fluid over a semi- 
infinite vertical plate has been studied when the free-strem velocity, surface 
mass transfer and wall temperature vary arbitrarily with time. The partial 
differential equations with three independent variables governing the flow 
have been solved using quasilinearization in combination with an implicit 
finite-difference scheme. The results indicate that the buoyancy para- 
meter, coupling parameter, mass transfer and unsteadiness in the free-stream 
velocity strongly affect the skin friction, microrotation gradient and heat 
transfer whereas the effect of microrotation parameter on the skin friction 
and heat transfer is rather weak, but microrotation gradient is strongly 
affected by it. The heat transfer is strongly dependent on the Prandtl 
number, the dissipation parameter and the variation of the wall temperature 
with time whereas the skin friction and microrotation gradient are weakly 
dependent on it. The buoyancy parameter causes an overshoot in the velocity 
profile. The magnitude of the velocity overshoot increases as the buoyancy 
parameter increases and it decreases as time increases. 


1. INTRODUCTION 


When the velocity of the fluid is small and the temperature difference between 
the surface and ambient fluid is large than the buoyancy effects on forced convective 
heat-transfer become important. The combined effect of forced and free convection 
over a heated vertical plate for Newtonian fluid has been studied by several investi- 
gators'-*, The flow and heat transfer behaviour of Polymeric fluids, colloidal fluids, 
real fluid with suspensions, liquid crystals and animal blood cannot be explained on 
the basis of Newtonian and non-Newtonian fluid theory. The theory of micropolar 
and thermomicropolar fluids was introduced by Eringen5~’, 
micropolar theory is given by Ariman et al,-°, 
the steady forced or free convection bound 
Recently, Jena and Mathur'® have studied t 
boundary-layer flow of a micropolar fluid from 
effects. It may be noted that the unsteady mixe 
in a micropolar fluid has not been studied so far 


An excellent review of the 
Several investigators'”!4 have studied 
ary-layer flow for micropolar fluids. 
he steady mixed convection laminar 
a vertical plate without dissipation 
d convection flow over a vertical plate 
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We have investigated the unsteady laminar incompressible boundary-layer flow 
of a micropolar fluid over a vertical flat plate when the free-stream velocity, mass 
transfer and the wall temperature vary arbitrarily with time. The effects of the surface 
mass transfer which varies arbitrarily with time, viscous dissipation and the Prandtl 
number have also been taken into account. The partial differential equations with 
three independent variables governing the flow have been solved numerically using a 
quasilinear finite-difference scheme. The results have been compared with Oosthuizen 
and Hart’, Gryzagoridis* and Jena and Mathur?’. 


2. GOVERNING EQUATIONS 


We consider the unsteady laminar incompressible boundary-layer flow of a 
thermomicropolar fluid past a vertical plate under the combined effect of forced and 
free convection. It has been assumed that the free-stream fluid temperature remains 
constant and the free-stream velocity, surface mass transfer and the wall temperature 
vary with time. Under the foregoing assumptions, the equations governing the flow 
can be written as°’””?°"?° 


Ux + vy = 0 ..-(1a) 
u, + ux + vuy = (Ue) + [(u + Ki))/P] uy + (Ay/P) Ny + 88 (T — Too) 
= (iD) 

N, + uNx + vNy = (y/P7) Nyy — (ki]Pj) [2N + wy] <es( 1c) 
T: + uTx + vTy = Pr * (u/P) Tyy + (1/Pcp) [7..Ny —TyNx] 

+ (Hep) | (aw + ky!2) wh + 2ky (N + wD 

+ ¥ Ny | : aid) 

The relevant initial and boundary conditions are 

u (x,y, 0) = us (x,y), ¥ (% Y, 0) = 1 Y) \ (2a) 
N (x, ys 0) = Mi (x, y), 74%; Vs 0) — Ti (x, y) J 
u (x, 0, t) = 0, v (x, 0, t) = vy (t) \ 
N (x, 0, t) = 0, T (x, 0, t) = TAL) r ...(2b) 


| 
u Lx, oo, t) = Ue (x, i: N (x, oo, t) = 0, Tix co, t)=T... J 


Here x and y are the distances along and perpendicular to the surface, Tespectts 
vely; ¢ is the time; u and v are the velocity components in the x and y cena 
respectively; N is the component of microrotation whose direction of patie ; 
the x — y plane; g is the acceleration due to gravity; Pp a nd 7 are ee: aa an 
temperature of the fluid; »,k: and y are the viscosity, vortex viscosity and Spl 
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gradient viscosity, respectively; 8 is the coefficient of thermal expansion; 1 is the 
micropolar heat conduction coefficient; cp is the specific heat of the fluid at constant 
pressure; j is the micro-inertia density; the subscripts t, x and y denote derivatives 
with respect to t, x and y, respectively; and the subscripts e and w denote conditions 
at the edge of the boundary layer and on the surface, respectively. The subscript i 
denotes values at the initial time t = 0 and 7, 1s a constant. 


It may be remarked that we have assumed that the microrotation N is equal to 
zero on thej boundary. The justification for using such a boundary condition is given in 
detail by Kirwan Jr’®. Here the free-stream velocity ue which vary with time can be 
expressed in the form 


Ue = Ugg 9 (t*), 1* = Uco tL. ay 


9 is an arbitrary function of the time /* representing the nature of unsteadiness in the 
external stream and has a continuous first derivative for t* > 0. 


On applying the transformations 


= (Ugg f2vEL)'l? y, EF = ¥, X = IL 

b = (2g & L)? f (E, n, t*) @ (t*) Ata 

u=u, F,f' = Fi, ue = tg 9 (t*) } 

V = —(Wo /2EL)!? o[f + 2 fe — nF] | 

N = [us 2E/(vL)p? s , ...(4b) 

Bice Too )(Two — od) =G 

(Tv — Too )(Tno — Too) = G 91 (t*) } 

f= J Fan + fof = A BMlg | 
... (4c 

A = —(Yy/Uo) (Rez/2)!/?, Rez = weg Liv } 


to (1), we find that (1a) is satisfied identically and (1b—d) reduce to 


(1 + M) F" + @ fF’ + 2 [9 g* (1 — F) — Fit] + 2€ 9 [Nys’ +AG] 


= 26 o [FFs — fe F’] ...(Sa) 
N3s" + N2o[fs' — sF] — N, [4 5 + p F’] — 2& N, s;* 
= 2 9 N, [Fee — fe 5’] ...(5b) 


Prt G" + 9 fG' +2€a [s Gg — 5G] — @ sG’ + Br[l + NJ2)¢2 F? 
+ 2M, Br [28 s + 9 F’/2? + 2 NsBrs'? — 2€ G,* 
= 26 9 [FG — fz G’] ..(5e) 


a 
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where 

N, = kyu, Ne = (voo/L) (P Jj] p) 

N3 = (Uo /L) (yP/p?), « = («,/cp) (Uco /L) 


A = Gr/Re, , Gr = g8 (Two — Too) L? p?|u? ...(6) 


e+ = ae 


Br= th, hcp (Two — T oo) ]- 


Here & and » are transformed coordinates; t* is the dimensionless time; u.. is 
the free-stream velocity at ** = 0; Lis the length of the plate; v is the kinematic 
viscosity; / and f are the dimensional and dimensionless stream functions, respectively; 
F(f'), s and G are the dimeasionless velocity, microrotation and temperature, respecti- 
vely. The parameters N,, N., Nz, « and Br are the coupling parameter, micro-intertia 
density parameter, microrotation parameter, micropolar heat conduction parameter 
and dissipation parameter, respectively; A is the buoyancy parameter; Gr is the 
Grashof number; and Pr is the Prandtl number. ff, is the surface mass transfer para- 
meter. If the normal velocity at the wall v, is selected in such a manner that (Vy/too) 
(Rez/2)'/? is a constant then the parameter A will be a constant. Hence the mass 
transfer parameter /, will vary according to (4c). A =] 0 according to whether there is 
a suction or injection. The subscripts € and t* denote derivatives with respect to § 
and t*, respectively. The prime denotes derivatives with respect to %. 


The transformed boundary conditions are given by 
F=0,s5s=—0,C = t*) aty = 0 
é jai ha \ for 1* > 0. Fae 


F>1l1,s-~0,G—>0 as 7) > co 


We assume that the flow is initially steady and then becomes unsteady for ¢* > 0. 
Hence the initial conditions for F, s and G at t* = 0 are given by steady flow equa- 
tions obtained by putting 


e(t*)=9(*) = 1, 0,= Fe=Se = GS, = 9 ---(8) 
‘ 


t t ‘ 
in (5) and also in (7) (for boundary conditions) and they are 
(1+ .N,) F’ + fF’ + 28 (N, 5’ + AG] = 26[ FFs — feF ...(9a) 
Nas” + N,[fs’ — sF] — N, [48s + F’] = 26N2 [Fse —fes'] _ ...(9b) 
Pr-) G’ + fG' + 2&x [s'Ge — sgG’'] — a sG’ + Br{t + N,/2] F” 


+ 2N, Br [2s + F]2]}? + 26 Ns Brs'? = 2&[FG~_ —f<G’]. 
...(9c) 


It may be remarked that the steady-state equations (9 (a)—(c)) reduce to that of 
Jena and Mathur’ if we replace 2& Nis’ by N,s’ in equation 9(a), —N.sF by N.sF and 
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—N, [4&s + F’] by —2EN, [2s + F’]in equation 9(b) and put « = Br = Oin equation 
9(c). The equations 9(a)—(c) also reduce to those of mixed convection for Newtonian 
fluids which have been studied by Oosthuizen and Hart* and Gryzagoridis‘ if we put 
N, = 0 and replace 2A G by AG in 9a), consequently, s = 0 and equation 9(b) 
becomes superfluous. 


The skin-friction coefficient at the wall is given by 
Cy = 2xy/[P (ue), J=(2Re "2 (1+ Me F. ...(10a) 
tf =O 


where 
Tw = [Cy + k,) uy + k, N]yn0- ...(10b) 


The heat-transfer coefficient in terms of Nusselt number is given by 
Nu = 2xqw/[ke (Two — Too)] = (2 Rex)'!? G ...(lla) 


where 
qu = [ke Ty + Be Nelyno. .».(11b) 


The couple stress coefficient is expressed in the form 
M = m,/[P (ue), x] = (Rex)-? Ns s* .e(12a) 
t 9 
where 


My = 7 (Ny)yeo. ...(12b) 


Here Cf, Nu and M are, respectively, the skin-friction coefficient, Nusselt number and 
couple stress coefficient; t», gw and my are, respectively, the shear stress, heat-transfer 
rate and couple stress at the wall; Bc: is the heat conduction parameter and k, is the 
thermal conductivity. 


3. RESULTS AND Discussion 


The equations (Sa)—(Sc) under boundary conditions (7) and initial conditions (9) 
have been solved numerically using an implicit finite-difference scheme with a quasi- 
linearization technique. Since the detailed description of the method is given in 
Bellman and Kalaba’? and Inouye and Tate’, its description is not repeated here. 
Computations have been a carried out for various values of the parameters 
A(-— 0.25 <A <€ 20), MOS5< N, & 13.5), Ns (0.5 << N, < 4.5), A(-O0.5<4 
< 0.5), Pr (0.7 < Pr<7.0) and Br (-—0.2 < Brg 0.2) with No=1.0, a =1.0,¢ = 0:2. 
¢, = 0.1, 4, = 0.1 and w* = 5.6, The unsteady free-stream velocity and wall 
temperature distributions considered here are given by 


PY = 1 +t, 9 (1) = 1 + 4 sin® (o* #*), 9 *) = 14 «78 
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where e, e, and e, are constants and «* is the frequency parameter. The effect of step 
sizes in y, & and f* directions and the edge of the boundary layer » 9 on the solution 
have been studied with a view to optimize them. Finally, the computations were 
carried out with An = 0.02, AE = 0.05, At* = 0.1 and no has been taken between 
4 and 8 depending upon the values of the parameters. The results presented here 
found to be independent of the step size An, AE, At* and the edge of the boundary 
layer yo. at least up to 4th decimal place. 


In order to asses the accurancy of method, we have compared our Nusselt 
number results for Newtonian fluids (V, = 0) with those of Oosthuizen and Hart® and 
Gryzagoridis*. We have also compared our skin-friction and heat-transfer results for 
micropolar fluids (Ni > 0) for steady flow (t* = 0) with those of Jena and Mathur’®. 


In both the cases, the results are found to be in good agreement and the comparison 
is shown in Figs. | and 2. 


0.8 


-/2 


x 


Nu Re 


° 2:0 4.0 6.0 8.0 10.0 
r 


=—1/3 


Fic. 1. Comparison of heat-transfer coefficient Nu Re. for 9 (t*) = 1.0, 91 (t*) = 1. 0, N, 


= N, = N; =3=A=Br=0. —-——~—, present method; 0, Oosthuizen and Hart; 4, 
Gryzagoridis. 


The results for the case 9 (t*) = 1 +e t™, « > 0 (accelerating flow) are given 
in Figs. 3-6 and those for the case 9 ((*) = 1 4+ «, sin? (* f*) (fluctating flow) in 
Fig. 7. 


The effect of buoyancy parameter A, mass transfer parameter A, Prandtl number 
Pr, coupling parameter \,, microrotation parameter Ns, variation of wall temperature 
with time o: (t*), distance § and dissipation parameter Br on the skin-friction, micro- 


rotation-gradient and heat-transfer parameters (Fi,—s,,—G,) are shown in 
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‘ ‘ a 
Fic. 2. Comparison of skin-friction parameter fel and heat-transfer parameter — Gi. for > (t*) 


=1.0, 9; (¢*) = 1.0, Ny = 0.1, N3 = 0.02, Pr = 9.0, A = 0.0,« = 0.0, Br = 0.0. ———, 
= 0; ———, 4 = 4; @, Jena and Mathur. 


Figs. 3-5. Figures 6(a)—(c) depict the effect of buoyancy parameter A on velocity, 
microrotation and temperature profiles. 


The results indicate that the skin-friction parameter FY, increases with time t* for 


A < 0 and decreases with it for A > 0 whereas the microrotation-gradient parameter 


—s,, increases with ¢* for values of A [see Figs. 3(a)-(b)]. For A < 0, the heat-transfer 
parameter —G’, increases with ¢* but for A > 0 it increases only after certain value of 
t* (Fig. 3(c)]. It has also been observed that for all r*, the parameters F’, —s‘' and 


—Gj, increase as A increases. Similar behaviour has also been observed by Jena and 


Mathur’ for steady case (** = 0). This is due to the fact that the buoyancy force 
(A > 0) gives rise to favourable pressure gradient which accelerates the fluid in the 
boundary layer and thereby increases the skin-friction, microrotation-gradient and 


heat-transfer parameters. The skin-friction, microrotation-gradient and heat-transfer 


parameters (F’, —s‘ , —G‘,) are reduced due to injection (A < 0) whatever may be 
the value of t* and the effect of suction (A > 0) is just the opposite. 

injection increases the momentum, microrotation and thermal bounda 
nesses which cause deceleration in the fluid and suctio 
3(a)-(b) also show that the 


This is because 
ry layer thick- 
n does the opposite. Figures 
skin-friction parameter F’ and microrotation-gradient 
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2.0 


0.5 





0 1.0 20 30 





Fic. 3. (b) microrotation-gradient parameter—s,, 


=Gy, 


eng 10 2.0 3.0 


Fia. 3 (c) Heat-transfer parameter =-G. for 9 (t*) =1 + et*, 9, (¢*) = 1.0,M = 1.5, Ns = 1.5, 


_—-— 4=-05, Pr=0.7; ———.———, A = 0.0, Pr = 0.7; 


= Q. = 0.5. 
Br = 0.0, & 6) A200, 


eee Ae OS err ONO A = 0.0; Pr=3 
Pr = 7.0. 
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parameter —s‘, decrease as Prandtl number Pr increases. Similar behaviour has also been 


observed by Wilks®. On the other hand, the heat-transfer parameter —G,, is found to 


increase with Pr (see Fig. 3(c)), because a large Prandtl number results in a thinner 
boundary layer with a corresponding large temperature at the wall and hence a large 
surface heat transfer. It has also been observed that the effect of Pr on the parameters 


F; and —s‘, is less as compared to the parameter — CH, 


Figures 4(a)-(c) show that any time /*, the skin-friction, microrotation-gradient 


and heat-transfer parameters (F,’, s’ 


w? 


—G' ) decrease as the coupling parameter N, 


increases. Similar behaviour has also been observed by other investigators'°?**°. The 
cause of this reduction is due to the thickening of momentum, microrotation and 
thermal boundary layers due to in the parameter N, which in turn decelerates the fluid 


in the boundary layer. As N; increases, the parameters F,, —s‘, and —G, decrease 
whatever may be the value of t*. It has also been observed from Figs. 4(a)-(c) that 


for § = 0, the parameters F’ , —s' and —G’ increase as time ¢* increases from 0 to 
3.0 but for § > 0, FY, decreases with t* and —s’ and —G" increase with it. We have 


also observed that as & increases, the parameters F’’ and —G' increase for A > 0 and 
they decrease for A <0. Similar trend has been observed by Jena and Mathur’® also. 
It is also clear from these figures that the parameters Fy, —s, and —G‘ increase with 
the variation of the wall temperature with t*. However, the heat-transfer parameter 


—G,; strongly affected by the variation of the wall temperature with t* whereas its effect 


on F’ and —s’ is rather weak. 


The effect of dissipation parameter Br on the parameters Foe ag and. — G sie 
depicted in Fig. 5. This figure show that the parameters FY and —s‘ increase as Br 
increases whereas the parameter —G‘, decreases with it. This behaviour is independent 
of the value of t*. For all values of Br, F;; decreases with ¢* (after certain t*) whereas 


~s,, increases with it. As 1* increases, the parameter — Gj, increases for Br < 0 and 
decreases for Br > 0, It has also beeen observed that the effect of Br is more pro- 
nounced on the parameter ~G\ than on the parameters F’ and —-s‘. 

w 


The effect of A on the velocity, 


microrotation and temperature is shown in Figs. 
6(a)-(c). Figure 6(a) shows that there 


is a velocity overshoot in F for buoyancy assisted 


~~ eS i—car,L”hmhCUhMTrhUrwemhC SCCLCUrCTLhLC™mCCUCrlUC™— ZO ee ae 
—— a - =r “es: =) e ) oo 
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“0 10 2.0 3.0 


e* 
Fic. 4. (c) Heattransfer parameter —Gi for 9 (t*) = 1 + et*°, Pr =0.7, A=0.0, A =5.0, Br=0.0. 


———,£ = 0.0, Ns = 1.5, 91 (t*) 1.0; ——————, §=0.5, Ns = 1.5. 91 (¢*) 
= 1.0; ———, & = 1.0, Ns = 1.5, o: (t*) = 1.0, --—O——-,, § = 0.5, Na = 4.5, 
g (t*) = 1.0, —-—A—— -, § = 0.5, Ns = 1.5, 91 (0%) = 1 + €2 0”. 
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3.0 
0 1.0 2.0 


Fic. 5. Skin-friction parameter Be , microrotation-gradient parameter mt ie and _heat-transfer 


parameter —G_, for »(t*)= 1+ et*?, 9, (t*)=1.0, N,; = 1.5, Ns = 1.5, A = 0.0, 


A= 5.0, Pr = 0.7, € = 0.5. ————, Br = 0.0, ———.———, Br = —0.2; ————, 
Br. = 0.2: 


flow (A > 0), and the velocity overshoot increses as A increases. However it decreases 
as t* increases. There is no velocity overshoot either for purely forced flow (A = 0) 
or buoyancy opposed flow (A < 0). The velocity overshoot is because buoyancy force 
(A > 0) gives rise to a favourable pressure gradient resulting in velocity which adds to 
the forced convection velocity. The buoyancy opposed flow (A < 0) gives to adverse 


ee —— ate 
—- . 
— . 
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Fic. 6. (a) Velocity profile in the x direction F. 
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pressure gradient which reduces the forced convection velocity. As the combined 


effect of buoyancy force A > 0 and forced convection force decrease with time which 
results in reduction in velocity overshoot with time. 


The microrotation profile —s and the temperature profile G are shown in Figs. 


6(b)-(c). It is observed that the profiles —s and G are significantly affected by the 
parameter A and the effect becomes more pronounced as time t* increases. 









The 
profiles —s and G become more steep as ¢* and A increase. 
.0 
ré) 1.0 2.0 3 5 ie 
Fic. 6. (c) Temperature profile G for (¢*) = 1 + et*?, 91 ") = 1:0; N, ae ne io 
Pr = 0.7, A = 0.0, Br = sce.) rs 0.0; ———. = 0: 


em, £* = 2,0. 
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The skin-friction, microrotation-gradient and heat-transfer parameters (F,,, —s,, 


i mie tg in 
—G' ) for oscillatory free-stream velocity 9 (t*) = 1 + «, sin? (#* ¢*) are shown 





0 1.0 ” 2.0 3.0 


Fic. 7. Skin-friction parameter Fo » microrotation-gradient parameter 5 and _heat-transfer 


parameter a for 9(t*) = 1 + e; sin? (w* ¢*), 9) (t*) = 1.0, N; = 1.5, Ng = 1.5, 
Pr= 0.7, A= 0.0, » = 5.0, Br = 0.0. ————, § = 0.5; ———-— § = 100. 


Fig. 7. It is clear from this figure that the parameters F\’, —s‘ and —G‘ oscillate 


as time f* increases but the oscillations are more for large &. 


4. CONCLUSIONS 
The skin friction, microrotation 


on the buoyancy parameter, coupling 
free-stream velocity. 


gradient and heat transfer are strongly dependent 
parameter, mass transfer and unsteadiness in 
The effect of microrotation parameter on microrotation gradient 


and heat transfer is comparatively 


» dissipation parameter and the variation of the wall 
temperature with time affect the heat transfer sig 


microrotation gradient are weakly affected by it. Buoyancy parameter induces over- 


eS ae 
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shoot in the velocity profiles. The magnitude of the velocity overshoot increases as the 
buoyancy parameter increases and it decreases as time increases. 
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The propagation of elastic disturbance ina nonhomogeneous medium due 
to radial pressure at the boundary of a spherical cavity has been studied. 
The problem admits of a closed form solution in some special cases only. 
The cases where no closed form solution is available have been investigated 
following Kromm’s technique. The solution has been obtained in terms of 
integral equation. Distribution of stresses have been shown graphically. 


1. INTRODUCTION 


The problem of the propagation of elastic disturbance in an infinite medium 
caused by a radial shock is the subject of interest of many investigators for its impor- 
tance in many practical fields viz. in seismology, in geophysical prospecting and in 
many other fields as has been discussed in detail by Neuber and Hahn}. We may 
think of a situation where there is spherical hole ina very large elastic medium and 
an explosion within the hole leads to release of huge amount of gas pressing the 
spherical boundary of the medium radially. Our problem is to study the propagation 
of elastic disturbance due to this radial pressure. The problem of this type for a 
homogeneous isotropic medium was apparently solved first by Jeffreys?. But as we 
know that all elastic media cannot be considered homogeneous and as it has been 
experimentally verified that there are plenty of media in which the elastic coefficients 
are functions of position, it is our motivation to study the problem for a more realistic 
medium, viz. a medium where the elastic co-efficients are not merely constants but are 
functions of position. The problem fora nonhomogeneous medium had also been 


considered by a number of investigators*-> but in all those Papers some restrictions 
were imposed to get the solution. 


In our present investigation we shall study the problem similar to those discussed 
above but with a modification that the boundary of the Cavity is disturbed suddenly 
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and thereafter maintained steadily, the cavity being free fromany other type of 
loading. With such a modification it is interesting to note that it has been possible to 
obtain complete solution of the problem for all values of the parameters involved. In 
our method of solution we have used the techniques followed by Kromm® and Goodier 
and Jahsman’. The expressions for stresses and displacement have been presented in 
the form of integral equations, which may be solved numerically. The expressions for 
the velocity of the wave front and also the values of the stresses and displacement at 
the wave front have been given. Finally the stress distributions for some particular 
cases have been shown graphically. 


2. FORMULATION OF THE PROBLEM 


Let the centre of the spherical cavity be taken as the pole of the spherical polar 
co-ordinate system (r, 6,9). In view of symmetry, the displacement components 
(u;, Ue, Ue), may be written as 


Up = u(r, t), Ue = Ue = O. 6 9 


The stress components *,r, tre etc. are connected with the displacement components 
(1) as, 


E Ou u 
tr = (+e) @ — 20) —sy[ 0-9) apes = | 


E fee) epi u 
Too = Toe = Ronis ar ae = | caek a 


99 = Tre = Tr, = 0,0 <  < 3. 


In (2) E is the variable Young’s modulus and othe constant Poisson’s ratio of the 
material. 


In view of (2) the equation of motion becomes 


0*u 
. Tr + 2 (trr — too) =P are: .-.(3) 


pP being the constant mass density. 


We assume that the Young’s modulus £ depends upon the position according 
to the law, 
E = E, (r/a)" (4) 
where parameter n may be any real number. 
Introducing nondimensional variables 


z = rlaand + = ct/a where c? = E,|P dest) 
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and taking into account of (2) and (4), the field equation (3) becomes, 
0*u 1 Ou 2no % ) ita ors AP Ou 6 
ee tt Dre + ( z2 oz” Ax? --«(6) 


in which 
oe. (1 + oc) (1 — 20) 
vty l—o f 
0 there is neither displacement nor velocity in the medium 


Assuming that at t 
the initial conditions of the problem may be put as, 
u(z,0)=0 
a 


2 u(z,0) = 0,1 < z < oc, 
As regards the boundary condition we assume that the boundary z = 1 is 


disturbed and thereafter the disturbance is steadily maintained. This may be expr- 
essed mathematically as, 
tr (1,7) = — PH(t),7 >0 ae FG 


where P is a constant and H (rt) is the Heaviside function. 
In addition to (8) the regularity condition requiring boundedness of the stresses 


and displacement for all values of z and t should be satisfied. 
Thus the solution of the problem reduces to the solution of (6) satisfying initial 
conditions (7), boundary condition (8) and also the regularity condition. 
3. SOLUTION OF THE PROBLEM 
In view of (7), Laplace transforms with respect to + of the field equation (6) and 


N | o> 
3| “to 
wD 
r 
= 


") 
| 


the boundary condition (8) yield 


d*ii 1 di 2no 
athe. fy Yate tee! ae 
Za tO) SF alee 2) 


and 


tr =—- —~— atz= ] 


where the bar over a function represents Laplace transform of the function with 


= 2andn > 2, 


respect to parameter p. 
Since the character of the solution of (9) depends on whether n < 2,n 

n > 2, we consider three cases, viz.,n < 2,n = . 

For n # 2, the general solution of (9) is, 


a = Az® K,, (pkz*) + B28 I, (pkz*) 
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where 


a = 1 — (n/2),B = — (n + 1))2,k = Afa 
and 








: 1 f(n + 1) 2 
m = | { ; vay ee ma|>0as0<c<4. 


l—o 


For n < 2, the argument of the Bessel functions in (11) is positive for Beto 
and tends to infinity as z > co, In this case the solution of (9) satisfying the regula- 
rity requirements is, 


ii = Az® Ky, (pkz*). ...(12) 
The arbitrary constant A = A (p) in (12) may be determined from (10). Finally we 
obtain, 


2Pa (1 + a) (1 — 20) 








i=— Paee Digest." Z(ntt)/2 K,, (pkz*) ra GY 
= Pz(n-3)/2 
rr = — —— [(1 — 0) (2 — 1) phz* Ki, (pkz*) 

+ {40 — (1 — oc) (n + 1)} Km (pkz*)] (14) 
| (n+3) /2 
00 = — [a (2 — n) pho* Ki, (pkz*) 

+ {2 — (n + 1) 6} Km (pkz*)] cis 


where 
D = (1 — s)(2 — n) pk K’ (pk) + {40 — (1 — a) (n + 1) }K,, (pk). 
Introducing the notations, 


a, =(1—)(2 —n) 
b, = 40 —(1 — oc) (n + 1) 


a, = a (2 = n) 
b= 2—(n+1)c 
@(z) = — ait (1 + 5) (1 — 20) z7(nti)/2 
Eo 
n—3 


—_— 


Y(z)= —Pz ? 
and using properties of Bessel function (cf. Abramowitz and Stegun*) we may rewrite 
(13), (14) and (15) as, 
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6 (z) Km (pkz*) 
i (z, p) = “Dp: A (p) ...(16) 
= v (z) Ai (Pp) 
rr (Z, P) = > a 2B (p) matias 
: v (z) a (p) 
oo (Z, P) = p A(p) .. (18) 
where 
A (p) = (a1 m + b,) K,, (pk) — ai Dk Kgs (pk) 
Ay (p) = (a1 m + bi) Km (pkz*) — a; pkz™ K,,41 (pkz*) (i = 1, 2). 
Let us set, 
9 (Pp) = (a,m+ b Da (Pk) — a, k Km4; (pk) ' 
| 
Hi (2, p) = 22) Kn (pk) | 
P | 
i 2) Km (pkz* ( 
B, (2, p) = 2) Sea, m + by) Rn (PREY a, ke Kina (pkz*)} | 
7 b (z ) Km (pkz* | 
ds (2, p) = Ita: m + b,) oo cntnsheee Fike (pkz*)}. | 
sateo) 


Then we may get (Erdelyi’) 
_ A(t —k) fam+b i a 
9 (0) = AS [a er + vera ie - co —veiain} 


— a, k (2 — YP (0, + YoFaTPH + (0, — vor} | 


hi) = vf 4G af — 4p [(o + vo® = Tye — (@ — Yor = oF 1 | dk 


bs (4) = 4 (2) paee eee 


0 
{( + Yot—1)" — (@ — Yo? — 1)m } 
— ak, (E2 — k2)0 {(o + Vato iyee 


+ (wo — vor Tm} | dé 
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bs (t) = (2) | a aed 


[+ Ve — w — Vora} 
= 1 ky B= PY (@ $ or = Tyme 


+ (@ — Vo? — ])mH \] dé .. (20) 


where 


©, = t/k, w = E/k, and k, = kz*. 


From (16), (17), (18) and (19) we have 


i (2, P)$ (p) = Vi (z, p) (21a) 
trr (z, P) 8 (p) = de (z, p) ...(21b) 
too (z, P)O(p) = bs (z, p). ...(21¢) 


Now taking inverse Laplace transform of (21a, b, c) we get, 


fu (z, 8) e(s — Ede =u (2, 2) (22a) 
Stn (z, E)@ (+ — &) dE = y (z, *) (2b) 
j reo (z, E) e(t — E) dE = Uy (z, 2). -..(22c) 


The integral equations (22a, b,c) may be reduced to comparatively simple forms 
suitable for numerical evaluation. Thus introducing the notations, 


gi (t) = ea (oN ag ag 1)? (0, — V/ OT De } 


_ 4 eed, fo, + fara I" 4(w0,— varaijme 


r-k, 


fi(Z,.5) = \ [(@ + ver= - (o— yor = 1" ] ak 


0 
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r 


i @, 1) = [ [Efe + yor 1m (o — vor — Hr} 


m 


— ak, (C2 + 20k)? 4 {( 4+ fo — 1)" 


+o vara Ye 


r- 


fs (2, t) = 


[seas 


2m 


{(o fav ja var iy} 


oe o> 


a a,k, (¢? + 20k,)-1/2 3 {(o a Vf ow? oe ])m+1 


+ (0 — yor Tym | lag 
(C= €& —k,) 


we obtain, 
g(t) = H(t — k) 9, (*) 
i (2, t) = 0 (z) H(t — ky) f, (z, 2) 
Ue (2, +) = b (z) H(t — ky) f, (z, 7) 
bs (2, t) = b (z) H(t — k,) fs (z, t) 
and eqns. (22a, b, c) take the form, 


j u(z,§)o(t1 +k — &)d— = H(t, — 8) 6 (2) 7, (2,41. K)> .5(23a) 


r 


i rr (2, &) 91 (% +k — A) dE = H (x, = he = +#) «..230) 


r 


free, Doi (ns +k — Od =H, — YO A( y+ W) (236) 


where 
t= 1t — kand 6 =k, —k. 


The integral equations (23b, c) are to be numerically solved to obtain the stresses 
Trr (2, t) and too (z, t). We shall follow the usual procedure of replacing the integral 
equation with approximately equivalent finite system of linear algebric equations, 


assuming that the unknown function is piecewise constant over the range of integra- 
tion. 


SL 
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Since the right-hand members of the equations (23b,c) contain H(t, — 41) as 
factor it follows that t,, (z,+,)=0 and ctoeo(z,t,)=0 for t, <8, and as 
8, = k(z* — 1), (k > 0, « > 0) it follows that 5 + co as z + co, This shows that 
the stress waves progress towards infinity in infinite time. If z and t momentarily, 
denote instantaneous radius of the wave front and the time at which this radius is 
attained, then the velocity of the propagation of the wave front is given by 


dz 1 
>= —_-— 


2 gals 
dt ae . 


Consequently velocity increases or decreases with z according asn > 0 or < 0. 


The values of the displacement and stresses at the wave front may be obtained 


as, 
u(z, 5; +) = oy ESD a 
e>0 
tre (z,5, +) = lim AGEAGS Tea: = b (z) 2°” 
e>0 
k ape 
too (Zz, 5, +) = lim eas = w (z) ee /2 
e>0 
where 
ete 
1= f 9,(6,+¢«+k — &) 4. 
§ 
1 
For n = 2, solution of (9) is given by, 
ec oe at Es ...(24) 
where 
m, m =4(-—3 + /9 + 4e,) 
and 


e, = A® p? + 2 — 4o/(1 — 9). er tS), 
We shall assume that the parameter p in Laplace transform is such that «, > 0 
in (25). So the appropriate solution will be # = B, 2™2. 


Applying the boundary condition (10) we can find B,. Finally the inversion 
formula of Laplace transform gives (Erdelyi’) 





= — 


Pa(1 + )(1 — 29) | s/s | g, (é) dé 
AE, (1 — 9) : 
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where 





BC) =f) —» F SOV =F), Ow) dy 


f(x) =e, r>y 
= 0 Oey 
y = logz(z > 1) 
n = (3 — 70)/2A (1 — a) 
v? = (17 — 330)/4A2 (1 — a) 
and J, is the Bessel function of the first kind and of order 1. 


The components of stresses are then computed as, 


Tre = — Pz? i 82 (5) dé 


ere Peet {gs (&) dé 


where 


82 (1) = 3 (sy) —v f B(V = Fy) J (w) dy 


> 


g(t) = g(t) —v j g (Vt? — »® — y) A (w) do 


ene tear (1 + 2c)(1 + 6) | 
OT ace ae A(L—efs, eee 








=-0,0<tT<yp 


3 (t) being the Dirac delta function. 


For n > 2, we find that « < 0 and therefore the argument of the Bessel func- 
tions is negative for p > 0, and then neither /,, nor Km is real valued. Moreover, the 
argument now tends to zero as zoo, Since at x = 0, J, (x) and K,, (x) have 
linearly independent singularities, the general solution (11) cannot satisfy regularity 
requirement. Thus we are compelled to abandon the regularity condition. This 
brings us to the difficulty of unique determination of two constants A and B in (11) 
with the help of single boundary condition (10), However, it has been pointed out by 
Sternberg and Chakravorty”® that for the possibility of a diverging wave ina physical 
domain B must be equal to zero and so the appropriate solution, in this case also, 
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will be given by (12). Thus proceeding as before we obtain the same expression for 
the displacement and stresses as in the casen < 2. 


4. NUMERICAL RESULTS 


In order to get some idea about the behaviour of the stresses with time at 


different positions of the medium for different values of the exponent we have 
computed the values of R = — "= and S = — > for different values of z, t+, and 
n corresponding to c = 1/3. In Fig. 1 and Fig. 3, the value of the nonhomogeneity 


0.6 ost 


0.5 0.5- 


04 





0.3 


R > 


‘0. 











0 2 4 6 a Lo 
q oe 
Fic. 1. Variation of t,- with time (for fixed 7). Fic. 2. Variation of t,; with time (for fixed 
position). 





Fic. 3. Variation of t gq with time (for fixed »). Fic. 4. Variation of teg with time (for fixed 
position). 
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parameter 7 has been kept fixed at1 and values of Rand Shave been evaluated 
respectively for different positions viz., z = 2,3 and4. In Fig. 2 and Fig. 4, the 
position z = 2 has been chosen and the values of Rand S have been calculated 
respectively for different values of, namely » = 0,0.5 and 1. The variations of 
stress components with position or with time is quite clear from the graphs. 
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